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SECTION  1 

BACKGROUND  AND  OUTLINE 


“The  bombs  bursting  in  air”  seldom  produce  craters,  but  most  other  explosions  do.  Hence,  anent 
weapons  effects  -  always  of  vital  concern  to  the  DoD  -  few  subjects  have  received  as  much  attention  as 
craters  and  cratering  (witness  CARES).  Two  decades  ago,  the  DoD  had  already  sponsored  at  least  25 
nuclear  explosive  events  and  150  of  conventional  type,  wholly  or  in  part  to  study  craters  (Cooper,  1976; 
a  partial  list),  plus  a  larger  number  of  “first  principles”  computer  simulations  of  such  events.  Indeed, 
as  recent  interest  in  tunnel  damage  illustrates,  cratering  remains  basic  among  weapon-effects  even  for 
conventional  weapons. 

Cratering  studies  (especially  theoretical  ones)  have  identified  several  processes  whereby  craters  form 
(“cratering  mechanisms”).  These  range  from  void  closure  (“compaction”),  the  creation  of  free-flying 
rock/soil  missiles  (ejecta  “throwout”  and  “fallback”),  and  permanent  displacement  during  shear  failure 
(“flow”) -which  occur  widely -to  more  subtle  and  specialized  processes  (e.g.,  “piping;”  at  wet,  atoll¬ 
like  sites,  rock  particles  in  slurries  are  forced  upward  through  fissures  by  geostatic  or  other  loads, 
whereupon  currents  in  otherwise  clear  water  carry  them  away).  Highly  specialized  mechanisms  abound, 
because  most  of  the  myriad  refinements  to  models  of  the  mechanical  properties  of  a  site  have  some 
effect  on  permanent  burst-induced  volume  changes  in  site  materials. 

Special  and  subtle  mechanisms  tend  to  be  more  speculative  and  less  effective  than  those  of  obvious 
widespread  occurrence  (such  as  compaction).  New  obvious  mechanisms  therefore  hold  relatively  great 
potential  for  explaining  crater  formation,  and  for  entering  decisively  into  cratering  theories.  Yet,  at 
this  late  date,  “new”  and  “obvious”  are  almost  contradictory  terms.  Both  apply  nonetheless  to  a 
nearly  universal  mechanism  (“obvious,”  after  the  fact)  that  surfaced  just  a  few  years  ago:  transport  of 
crater  volume  to  very  great  distances  by  elastic  waves  of  low  (even  seismic)  amplitude.  Accordingly, 
crater-volume  budgets  need  to  specify,  at  the  least  (i.e.,  ignoring  special  subtle  mechanisms),  the 
contributions  to  total  volume  from  throwout,  compaction,  flow,  and  elastic  volume  transport. 


The  contribution  made  to  crater  volume  by  elastic  volume  transport  is  the  main  subject  of  this  report. 
Further,  we  consider  throughout  only  craters  to  which  throwout  contributes  nothing,  a  simple  limit 
suggested  by  DNA’s  PEACE  exploration  of  Pacific  craters  (which  stimulated  this  work)(Wardlaw, 
1987). 

The  report’s  next  section  establishes  elastic  volume  transport  as  a  major  cratering  mechanism.  Then, 
drawing  as  needed  on  knowledge  of  compaction,  flow  and  elastic  deformation,  and  on  exact  fields  of 
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motion  in  which  one  or  more  of  these  three  processes  occur,  we  calculate  and  tabulate  fractions  of  total 
crater  volume  that  might  result  from  them  in  various  media.  The  last  fraction  (that  for  elastic 
transport)  runs  from  about  £  to  §,  but  our  results  are  approximate;  for  axisymmetric  fields  in  which  all 
three  mechanisms  operate,  the  equations  of  continuum  motion  have  so  far  defied  exact  solution. 
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SECTION  2 

TRANSPORT  OF  VOLUME  BY  ELASTIC  WAVES 
2.1  ELEMENTARY  FIELDS. 

We  recall  first  the  case  of  a  flat  elastic  halfspace  uniformly  loaded  by  a  pressure  pulse  of  finite 
duration.  The  well-known  result:  The  surface  receives  a  permanent  downward  displacement  equal  to 
the  impulse  delivered  per  unit  area,  divided  by  acoustic  impedance  (i.e.,  longitudinal  wavespeed  times 
density)  —  and  the  initial  and  final  densities  of  material  are  identical.  This  result  (no  density  change) 
comes  about  as  follows:  The  volume  through  which  the  surface  drops  is  equal  precisely  to  the  volume 
of  material  removed  by  compression  in  the  wave  traveling  downward  through  the  medium.  Thus,  on 
asking  where  the  volume  of  this  “crater”  lies,  we  find  it  moving  down  endlessly  at  constant  speed,  ever 
farther  from  the  “cratered”  (depressed)  surface. 


Turning  next  to  the  case  of  a  spherical  cavity  in  a  homogeneous  isotropic  elastic  wholespace,  with 
pressure  acting  uniformly  over  the  cavity  wall,  let  that  pressure  grow  to  some  new  fixed  value.  The 
cavity  radius  then  grows  too.  When  the  system  stabilizes  at  the  new  cavity  pressure,  however,  the 
density  throughout  the  medium  is  again  identical  to  its  initial  value.  As  in  the  case  of  the  flat 
halfspace,  a  load  impressed  on  the  medium’s  surface  produces  no  net  change  in  volume  strain;  once 
more,  a  cavity-volume  increase  is  made  possible  by  compression  in  a  wave  that  runs  away  from  the 
cavity  at  constant  speed,  forever.  Nor  does  similarity  to  the  halfspace  field  end  there:  Below  a  plane 
through  the  center  of  symmetry,  the  added  “crater”  volume  (half  the  change  in  cavity  volume)  is  again 
proportional  to  the  impulse  that  the  pressure  change  delivers  (to  the  medium)  through  that  plane.  The 
impulse  in  question  remains  finite,  be  it  noted,  even  though  the  pressure  change  lasts  for  all  time  at  the 
wall.  [Physically:  A  vessel  can  confine  fluid  during  and  after  a  pressure  buildup  because  tensile 
stresses  develop  in  vessel  material;  on  either  side  of  a  cross-sectional  plane,  the  tensile  force  in  the 
material  offsets  the  force  on  the  vessel  due  to  fluid  pressure,  with  the  sum  of  the  forces  becoming  zero 
at  equilibrium.  In  holding  cavity  pressure  forever,  the  wholespace  acts  like  a  vessel  with  very  thick 
walls:  In  planar  cross-section  through  the  cavity’s  center,  tractions  in  the  material  are  indeed  tensile 
(and  normal  to  the  plane).]  These  results,  elementary  but  seldom  brought  up,  are  derived  in 
Appendix  A. 

Linear  elasticity,  however,  does  not  approximate  well  the  behavior  of  material  near  the  craters  that 
explosions  produce.  Indeed,  it  does  not  allow  permanent  craters  to  form;  as  unloading  proceeds,  the 
ground  surface  regains  its  pre-shot  shape.  True,  in  a  wholespace,  permanent  pressure  in  the  cavity 
would  leave  it  permanently  enlarged,  but  (as  big-hole  decoupling  demonstrates)  by  a  small  fraction  of 
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the  volume  that  an  actual  explosion  produces.  Rather,  craterologists  have  invoked  the  inelastic 
processes  of  flow  and  compaction,  mainly,  to  account  both  for  permanent  cratering  and  the  growth  of 
cavities  in  deeply  buried  bursts.  Compaction  aside,  in  an  elastic-plastic  medium,  material  deforms  for 
a  while  at  the  maximum  shear  stress  it  can  sustain  (i.e.,  it  inhabits  states  of  shear  failure),  but  with 
mean  stresses  of  much  higher  amplitude;  it  deforms,  therefore,  almost  in  fluid  (i.e.,  hydrostatic)  states 
of  stress -it  flows.  Even  in  a  wholespace,  however,  stress  levels  decay  with  distance  from  shot-point 
(and,  near  the  cavity,  with  time;  cavity  pressure  drops  as  it  expands);  hence,  shear  failure  does  not 
occur  beyond  some  finite  distance  (the  “elastic-plastic”  radius).  Outward  motion  of  the  spherical  shell 
of  material  at  that  radius  increases  the  volume  contained  within  it,  and  hence  adds  to  the  cavity’s 
volume.  Moreover,  since  deformation  is  elastic  at  all  larger  radii,  that  added  volume  again  resides  in 
material  compression  in  an  elastic  wave,  ever  more  remote,  that  moves  outward  at  constant  speed.  At 
any  given  distance  of  propagation,  however,  the  amplitude  of  the  wave  is  enhanced  (often  greatly)  by 
the  incidence  of  shear  failure  within  the  elastic-plastic  radius -one  major  aspect  (among  several  that 
we  discuss)  in  which  cavity  growth  and  crater  growth  are  closely  related.  Next,  therefore,  we  recall  in 
more  detail  what  goes  on  near  the  cavity,  and  add  compaction  to  the  combined  effects  of  flow  and 
elasticity. 

2.2  PRINCIPAL  STRESS,  SHEAR  STRESS,  AND  SPHERICAL  MOTION. 

The  detonation  front  in  a  conventional  explosion  is,  ideally,  a  sharp  shock,  and  can  be  treated  as  such 
for  most  weapons  effects  purposes  (by  far);  the  same  is  then  true  (almost  always)  of  the  wavefront 
driven  outward  in  adjacent  material.  In  a  nuclear  burst,  radiation-hydrodynamic  interaction  may 
occur  at  first  between  explosive  and  medium,  but  a  sharp  shock  nonetheless  develops  in  adjacent 

material,  and  soon  (on  the  scale  of  crater  or  cavity  size).  In  isotropic  material,  the  shock  causes 

straining  only  along  its  direction  of  travel  (“uniaxial  strain”).  Strain  in  that  direction  is  a  principal 

strain,  as  are  the  null  “lateral”  strains  in  directions  orthogonal  to  the  shock.  The  corresponding 

components  of  stress  (normal  stresses)  are  also  principal  values,  but  (with  rare  exceptions)  lateral 
stresses  do  not  equal  zero.  In  fact,  all  measurements  say  that  lateral  stresses  can  be  as  large  as  axial 
stress  (fluids),  but  are  smaller  in  solids;  i.e.,  the  pressure  applied  to  the  ends  of  a  solid  cylinder  to 
compress  it  axially  exceeds  the  pressure  that  must  then  be  applied  to  keep  it  from  expanding  radially. 

The  directions  in  which  principal  stresses  act  -  principal  stress  “axes”  -  are  special  m  that,  on  a  plane 
normal  to  any  given  axis,  shear  stresses  vanish;  across  that  plane,  the  force  on  material  (per  unit  area: 
“traction”)  is  axis-directed  (i.e.,  normal  to  the  plane).  Across  any  other  plane  (neither  parallel  nor 
orthogonal  to  the  shock  direction,  in  the  present  instance)  tangential  components  of  traction  (shear 
stresses)  are  linear  combinations  of  differences  between  principal  stresses.  Hence,  if  all  principal-stress 
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differences  vanish  (principal  stresses  equal,  as  in  a  fluid),  then  so  do  shear  stresses  (as  in  a  fluid)  —  and 
the  reverse  also  holds.  Further,  since  there  is  a  finite  limit  to  the  size  of  the  shear  stress  that  can 
develop  in  any  real  material,  its  principal-stress  differences  are  also  limited  in  size  (at  the  limit  or  not, 
half  of  the  largest  difference  between  a  pair  of  principal  stresses  equals  the  largest  shear  traction  found 
on  all  planes  through  a  given  spatial  point).  Less  formally  put,  matter  does  not  have  the  “strength”  to 
develop  larger  shear  stresses  or  principal-stress  differences  than  this  material-specific  limit  allows; 
driven  to  its  limit,  it  deforms  (“flows”)  instead,  without  limit. 


Principal  stresses,  being  normal  components  of  traction,  can  change  during  shear  failure,  but  with  the 
largest  principal-stress  difference  fixed.  Hence,  in  spherical  fields,  two  of  whose  principal  stresses  are 
always  equal,  a  change  in  any  one  of  them  must  result  in  the  same  change  in  the  other  two  during 
shear  failure  (the  mean  of  the  squares  of  their  differences  must  not  change).  With  shear  stress  at  its 
limit,  changes  in  stress  are  thus  isotropic,  precisely  as  they  are  in  a  fluid.  Above  the  level  that  suffices 
to  generate  that  limit  at  the  cavity  wall,  pressure  on  the  wall  therefore  tends  (as  in  a  fluid)  to  be 
imparted  to  the  entire  region  of  limiting  stress,  or  “shear-failure”  region  (a  region  enlarged  by  the 
field’s  dynamics).  Intuitively,  such  an  outward  projection  of  wall  pressure  should  increase  the 
displacement  of  surrounding  material,  adding  to  cavity  growth.  Strong  support  is  lent  to  that  notion 
simply  by  examining  the  equation  (from  mass  and  momentum  conservation;  Hopkins,  1960)  that 
governs  spherical- fields-in-general  (whether  elastic  or  not),  namely 

P\  =  -  [d(Tr/dr  +  (2/r)(err  -  <r0)]  (2.1) 

where  p  and  ar  are  the  density  and  radial  acceleration,  respectively,  of  material  at  a  point  of  radius  r 
from  the  center  of  spherical  symmetry;  <rr  and  a6  are  the  radial  and  hoop  components  of  stress  at  that 
point,  respectively  (each  being  a  principal  stress,  and  positive  in  compression ). 

Specifically,  in  the  outgoing  shock  from  an  explosion,  stress  in  the  direction  of  the  shock  (a  principal 
stress;  above)  becomes  radial  stress;  the  lateral  principal  stresses  —  hoop  stresses  —  are  equal.  Thus,  <rT 
starts  out  compressive  and  hence  positive,  and  <Tq  <  ar  (above),  so  that  ( crr  —  a $)  is  at  first  >  0. 
Plainly  [Equation  (2.1)],  the  larger  the  value  of  the  principal-stress  difference  {<rr-(Te\  the  smaller  the 
radial  acceleration  of  just-shocked  material.  A  limit  to  the  size  of  shear  stress,  which  now  becomes  a 
limit  on  T|crr  —  cr^  |  (above),  therefore  acts  to  increase  the  radial  acceleration  of  that  material.  In 
elastic  media,  dajdr  (like  wall  pressure  and  <rr)  is  also  positive  at  and  near  the  shock  ( ar  falls  on 
moving  inward  from  it;  Cooper,  1966);  by  enhancing  outward  transmission  of  wall  pressure,  shear 
failure  causes  dar/d r  to  decrease,  and  therefore  again  acts  to  increase  ar. 

Note  too  that  as  just-shocked  material  moves  outward,  its  hoop  strains  (both  equal)  become  more 
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tensile  (circles  of  particles  about  the  center  of  symmetry  expand).  Its  radial  strain,  however,  tends  to 
become  more  compressive;  if  a  thin  spherical  shell  of  material  were  to  move  outward  at  constant 
volume,  its  rate  of  change  of  radial  strain  would  be  compressive  and  twice  as  large  as  that  of  hoop 
strain.  Indeed,  for  radial  strain  er  to  be  expansive  at  a  greater  rate  than  hoop  strain  eg,  volume  strain 
would  have  to  grow  at  more  than  three  times  the  rate  of  hoop  strain.  Until  this  last  condition  is 
reached  where  'signifies  time-rate-of-change  in  a  material  element),  ar  in  a  real  material 

may  not  decrease  faster  than  ae.  In  isotropic  linear  elastic  media  (as  it  happens)  -&r  >-&6  (hence 
&r-6e<Q)  if  and  only  if  -er>-e0.  This  need  not  be  so  for  other  substances;  e.g.,  given  shock-initiated 
shear  failure,  -&r  can  at  least  equal  -&e  during  the  period  of  flow  in  which  -er<-e0.  Yet,  the  later 
stage  of  spherical  motion  in  which  -tr>-ee,  is  reached  sooner  without  shear  failure  overriding  linear 
elasticity.  That  is,  by  admitting  large  positive  values  of  <rr-<re  in  the  shock  state,  linear  elastic 
behavior  brakes  a  medium’s  outward  motion,  thus  increasing  ie  (making  it  less  tensile)  and  reducing 
er,  and  so  gives  rise  more  readily  to  the  condition  -er>-e6.  Reaching  it  (-er>-e0)  sooner  during 
elastic  deformation  than  during  flow  thus  goes  hand  in  hand  with  less  cavity  growth  in  the  elastic  case. 

Stated  in  other  (more  common)  terms:  For  a  wide  class  of  materials,  increments  in  shear  stress  are 
proportional  to  shear-strain  increments.  The  ratio  of  the  former  to  the  latter  defines  a  shear  modulus 
G  that  generally  depends  on  the  stress-strain  state  from  which  the  increment  occurs,  and  may  depend 
on  an  element’s  entire  stress-strain  history  in  reaching  that  state  (although  for  linear  elastic  materials, 
G  is  constant).  G  measures  the  material’s  resistance  to  shear  deformation;  the  larger  the  value  of  G, 
the  greater  the  shear-stress  increment  required  to  produce  a  given  increment  of  shear-strain.  And, 
during  shear  failure,  G  can  be  zero  (Appendix  B).  With  the  shock  at  a  given  radius,  however,  a  larger 
cavity  implies  larger  values  of  er  and  -e0  (on  average  over  disturbed  material),  and  hence  a  larger 
shear  strain  [of  magnitude  \  \  er  -  ee  |  in  spherical  fields]  -  an  increase  unopposed  by  resistance  to  more 
shear  when  G=0,  and  resisted  more  and  more  strongly  for  larger  and  larger  G.  Small  values  of  G  go 
hand  in  hand  with  prodigious  cavity  growth. 

2.3  DISPLACEMENT  AROUND  SPHERICAL  CAVITIES:  LINEAR  ELASTICITY  AND 
QUASISTATIC  ELASTO-PLASTICITY. 

Linear  elastic  media  can  have  a  wide  range  of  G-values  even  if  they  have  the  same  density  p  and 
longitudinal  (radial)  wavespeed  cp  (hence  impedance  pcp  and  constrained  modulus  M=pcp).  They  can 
therefore  furnish  a  quantitative  sense  of  the  dependence  of  cavity  growth  on  G.  In  particular,  if  the 
cavity  pressure  jumps  suddenly  to  a  constant  level  P  beyond  overburden  (i.e.,  beyond  mean  stress  in 
undisturbed  material),  then,  after  motion  stops,  the  wall  will  have  been  displaced  outward,  relative  to 
the  cavity’s  radius  for  P=0,  by  a  fraction  |P/G  of  that  radius  [Appendix  A,  Equation  (A.ll)].  Thus, 
as  G-*0  [elastic  fluid;  Poisson’s  ratio  (i/)=|],  the  cavity  grows  without  limit  (a  correct  conclusion,  even 
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though  the  linearization  employed  in  reaching  it  is  acceptable  only  if  ip/G«l).  Moreover,  as  the 
field  of  motion  develops,  is  positive  in  just-shocked  material  until  the  shock  has  traveled 

-M/G-1  initial  cavity  radii.  At  the  cavity  wall,  which  oscillates  about  its  final  radius,  frr-&Q  remains 
>0  for  a  time-interval  determined  (in  units  of  the  time  taken  by  the  shock  to  travel  one  cavity  radius) 
by  G/M  alone,  and  hence  on'  v  alone;  that  interval  never  ends  when  (G=0),  but  decreases 

monotonically  to  2.03  units  of  time  at  i/=±  and  f  units  at  v=0  (G=|M).  Elastic  or  not,  the  shock- 
state  with  which  an  element’s  motion  begins  has  <rf  —  a q  >  0,  and  <xr  — 0^  starts  to  decrease  when  -<xr 
if  linear  elasticity  holds  up  to  a  given  “shock  stress”  (the  value  of  <j f  in  the  shock  state),  then 
in  just-shocked  material  approaches  zero  in  the  elastic-fluid  limit  [  i/(Poisson’s  ratio)=j],  and 
increases  steadily  to  <rr  itself  as  v  decreases  to  zero. 

Insight  is  gained  more  directly  into  the  effects  of  flow  on  cavity  size  by  considering  quasistatic  buildup 
of  cavity  pressure  in  a  homogeneous  isotropic  medium  whose  shear-stress  amplitude  cannot  exceed  a 
fixed  constant  limit  (von  Mises  limit),  but  is  otherwise  linear-elastic.  For  spherical  fields,  said  limit  on 
shear  stress  is  evidently  (above)  equivalent  to  the  constraint  |  <rr  |  <  Y,  where  the  strength  Y  is 
a  material-specific  constant.  In  this  case,  the  known  exact  field  (Hill,  1956,  p.  98-102)  makes  P/Y 
critical  to  determining  the  wall’s  outward  displacement.  For  one  thing,  shear  failure  occurs  around  the 
cavity  if  and  only  if  P>§Y.  Then  (for  example),  with  v~\  and  K/Y=50  (K  =  bulk  modulus^  M-§G), 
pressurizing  the  cavity  causes  its  radius  to  increase  by  a  factor  of  ^  if  P/Y=2.06,  and  a  factor  of  2  if 
P/Y=2.61;  with  K/Y=500  (and  i/= |),  cavity  radius  increases  by  factors  of  §  and  2  for  P/Y-values  of 
3.49  and  4.09,  respectively. 

Of  particular  interest  in  these  cases  (z/=^  and  K/Y-values  of  50  and  500)  is  the  source  of  the  eightfold 
volume  increase  that  occurs  when  the  radius  doubles  —  an  addition  to  initial  cavity  volume  V0  of 
7xV0:  With  K/Y=50,  material  fails  in  shear  out  to  5.25  initial  cavity  radii,  a  boundary  at  which  the 
outward  displacement  of  material  opens  up  a  volume  of  2.925  x  VQ;  the  same  outward  transmission  of 
cavity  pressure  (due  to  fluid-like  behavior  of  failing  material)  responsible  for  that  displacement  also 
compresses  failed  material,  shrinking  the  volume  it  occupies  by  4.075  x  V0.  For  K/Y=500,  the  volume 
(7  x  VQ)  added  to  the  initial  cavity  has  components  2.71  xV0  due  to  outward  displacement  of  material 
at  the  elastic-plastic  boundary  (initially  at  11.06  radii)  and  4.29  xV0  due  to  compression  of  failed 
material.  The  profound  effect  of  shear  failure  on  cavity  volume  can  be  seen  by  considering  buildups  of 
cavity  pressure  identical  to  these,  but  with  the  medium  strong  enough  (yet  otherwise  identical)  to 
preclude  shear  failure  at  those  pressures  (Y  >  §P).  In  the  resulting  purely  elastic  fields,  outward 
displacement  of  the  same  material  shells  (initially  at  5.25  and  11.06  cavity  radii)  open  up  volumes  that 
equal  just  .070  x  V0  (for  P/K=2.61/50)  and  .010  xVQ  (P/K=4.09/500).  Moreover,  except  for  nonlinear 
effects  (very  small  here)  of  finite  displacement,  volume  strains  are  zero  in  spherical-elastic  fields  around 
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cavities;  compression  of  surrounding  material  adds  nothing  to  cavity  volume.  And  thereby  hangs 
another  tale  of  central  importance  later  in  this  report:  As  already  noted  (Section  2.1),  the  medium’s 
stress  (hoop  stress,  is  tensile  over  a  plane  through  the  center  of  the  cavity  if  the  medium  is  elastic; 
when  shear  failure  occurs,  however,  is  compressive  over  the  inner  exp(-i)  (=.846)  of  the  distance  to 
the  elastic-plastic  boundary. 

Subsequent  quasistatic  removal  of  pressure  from  the  cavity  wall  (“unloading”)  raises  further  basic 
points.  First,  the  wall  will  surely  move  inward  as  the  pressure  on  it  drops.  Near  the  wall,  that  makes 
e6  (hence  <r0)  less  and  and  less  tensile  (more  compressive)  even  as  <rr  becomes  less  compressive,  so  that 
ar-cre  falls  from  the  value  Y  which,  prior  to  unloading,  obtains  over  the  shear-failure  region. 
Unloading  therefore  starts  as  an  elastic  process,  remaining  so  unless  and  until  it  proceeds  so  far  that 
<rr  -  ae  reaches  the  value  -Y  (shear  resistance  to  inward  displacement  then  does  not  increase  on  further 
unloading).  In  particular,  material  will  fail  in  the  latter  stages  of  unloading  if  P>|Y,  as  in  the  two 
examples  above  (u=\-,  cavity  radius  doubles;  P/Y=2.61,  4.09).  More  significant,  though,  is  the  small 
displacement  that  attends  purely  elastic  unloading,  which,  in  one  case  (P/Y=2.61),  continues  until  the 
pressure  has  dropped  at  the  wall  to  less  than  half  of  P.  At  that  point,  the  decrease  in  cavity  radius 
amounts  only  to  2§%  of  the  growth  that  took  place  during  flow  (when  shear-strain  increases  went 
unresisted).  For,  flow  stops  when  shearing  is  resisted  ( elastic  unloading);  material  “freezes.”  True, 
flow  will  occur  on  still  further  unloading,  but  that  fact  holds  only  minor  interest  here:  In  explosions, 
wall  pressure  drops  a  great  deal  as  the  cavity  expands;  eventually,  cavity  growth  stems  not  from 
pressure  on  the  wall  but  from  the  outward  velocity  already  imparted  to  material  around  the  cavity. 


2.4  SPHERICAL  CAVITIES:  DYNAMIC  EFFECTS  IN  ELASTO-PLASTIC  MEDIA. 

Indeed,  far  from  being  quasistatic,  cavity  formation  in  an  explosion  is  largely  a  matter  of  wave 
propagation;  as  the  process  unfolds  in  time,  it  presents  us  with  some  key  facts  and  phenomena  that  do 
not  enter  the  quasistatic  limit.  In  particular,  in  an  explosion,  the  cavity  reaches  static  equilibrium 
(apart  from  minor  damped  oscillations  about  its  final  state)  in  a  time  inversely  proportional  to  the 
medium’s  longitudinal  wavespeed;  dome-salt,  in  which  near-spherical  free-field  motion  has  probably 
been  best  measured,  has  a  wavespeed  of  4.2  km/sec,  and  cavities  driven  by  tamped  bursts  therein 
stabilize  within  4  ms/T1/3.  For  the  simple  elasto-plastic  material  above  (relevant  for  dome-salt), 
stabilization  comes  about  as  shear  failure  ceases  to  occur  in  material  around  the  cavity -the  result  of 
decay  of  both  cavity  pressure  and  shock  stress.  Material  there  unloads  elastically,  regaining  the 
capacity  to  resist  the  changes  in  shear  strain  whose  earlier  unopposed  development  (in  states  of  shear 
failure,  starting  with  the  initial  shock  state  itself)  made  large  outward  displacements  possible:  The 
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medium  stops  flowing  (failing  in  shear);  it  “freezes”  by  turning  elastic  as  stress  levels  drop,  just  as  it 
does  during  quasistatic  relief  of  cavity  pressure  (above). 

Dynamically,  of  course,  material  is  disturbed  at  that  time  only  out  to  the  shock’s  finite  radius;  on  that 
region,  radial  and  hoop  stresses  do  not  exhibit  quasistatic  profiles  —  except  that,  at  the  cavity  wall, 
radial  stress  must  equal  cavity  pressure  at  all  times  in  both  cases.  More  important,  an  appreciable 
fraction  of  the  work  done  on  the  medium  resides  in  its  kinetic  energy  (zero  quasistatically).  As  a 
result,  material  continues  to  move  outward  well  after  —  & r  first  exceeds  —  &q  (the  start  of  elastic 
unloading),  because  of  matter’s  inertial  resistance  to  changes  in  velocity.  Thus,  the  cavity  at  first 
expands  beyond  its  ultimate  (equilibrium)  size.  Then,  having  developed  a  finite  inward  velocity  in 
returning  to  that  size,  it  moves  inward  still  further  (inertia  again);  that  gives  rise  to  compressive  hoop 
strains  and  stresses  (relative  to  the  equilibrium  state)  in  near-cavity  material,  which  therefore  slows 
and  reverses  direction  -  etc.  There  thus  ensues  a  series  of  nearly-free  oscillations  of  the  cavity  wall 
( nearly ;  cavity  gases  affect  it  slightly),  diminishing  in  amplitude  because  the  waves  it  generates  in  the 
medium,  even  if  elastic,  carry  energy  away  from  the  cavity. 

The  oscillations  are  in  fact  elastic  unless  unloading  is  so  strong  in  the  first  of  them  that  aT  -  a0  drops 
all  the  way  from  Y  (in  the  states  of  shear-failure  that  start  with  the  shock)  to  — Y;  this  possibility  is 
explored  in  Appendix  B  (which  treats  the  subject  in  more  detail),  not  because  it  alters  much  the 
picture  drawn  here  of  cavity  stabilization,  but  because  it  explains  why  locked-in  stresses  (which  are 
compressive)  are  limited  in  amplitude  to  Y.  The  larger  fact  to  grasp  is  that  the  changes  in  stress  which 
cause  near-cavity  material  to  keep  reversing  direction  are  largely  (often  entirely)  elastic,  and  so  take 
place  with  far  less  radial  displacement  than  shear  failure  visited  on  the  material;  material  freezes  when 
shear-failure  stops  (which  happens  just  before  it  first  reverses  direction;  Appendix  B.4).  The  reason 
(again):  In  shear-failure  states,  the  growth  of  shear  strain  is  not  resisted;  in  elastic  states,  it  is. 

2.5  SPHERICAL  CAVITIES:  COMPACTION. 

Like  flow,  in  which  the  work  done  to  create  non-recover  able  (i.e.,  inelastic)  shear  strain  is  converted  to 
heat,  any  inelastic  decrease  in  volume  during  hydrostatic  compression  (“compaction”)  converts  work  to 
heat  (mainly  during  collapse  of  pores).  Likewise,  in  spherical  fields  during  flow,  reduced  resistance  to 
shearing  permits  greater-than-elastic  radial  compression  at  a  given  radial  stress  (Appendix  B.2),  and 
the  loose  packing  of  grains  or  blocks  of  solid  that  leads  to  compaction  (or  even  loose  molecule-packing 
prior  to  transitions  to  denser  phases)  also  enhances  radial  compression  at  that  stress  (cavity  pressure, 
not  medium- properties,  sets  radial  stress  at  first).  Either  way,  enhanced  radial  compression  hastens 
early  cavity  growth  (under  a  given  pressure,  the  softer  the  medium  the  faster  it  moves).  But-  whereas 


stress  is  transferred  outward  from  the  cavity  wall  through  material  that  flows  (above),  the  wall’s 
motion  is  projected  outward  through  material  that  compacts . 

The  reason  lies  in  the  relatively  large  increase  in  wavespeed  typical  of  compaction.  That  is,  compacted 
material  tends  to  be  much  stiffer  in  compression,  and  especially  in  subsequent  expansion,  than  in  its 
initial  loose  state.  Wavelets  that  signal  slowing  of  the  wall  can  therefore  traverse  disturbed  material 
while  the  outgoing  shock  (the  limit  of  disturbance)  moves  just  a  small  fraction  of  the  cavity’s  radius. 
At  the  same  time,  great  stiffness  in  compression/expansion  means  that  when  the  wall’s  radius  changes 
by  a  given  fraction,  the  volume  through  which  it  moves  is  equaled,  roughly,  by  the  volume  through 
which  any  spherical  sheet  of  compacted  material  moves.  That  is,  the  whole  region  of  compacted 
material  tends  to  move  as  a  spherical  piston  of  growing  mass:  As  the  wall  slows,  it  all  slows. 

By  the  same  token,  if  any  sheet  of  compacted  material  were  to  slow,  so  would  the  wall  -  and,  were  it 
not  for  shear  failure ,  that  would  actually  happen  before  the  cavity  could  grow  appreciably.  For, 
conversion  of  the  medium  to  stiff  elastic  material  via  compaction  would  entail  rapidly  increasing 
resistance  to  the  shear-strain  changes  implicit  in  outward  radial  motion  (Sections  2.2,  2.4;  Appendix 
B.3).  High  porosity  and  high  strength,  however,  are  strange  bedfellows.  Moreover,  material  does 
harden  (grow  stronger)  as  it  compacts,  but  only  toward  the  strength  it  has  when  consolidated.  Thus, 
accounting  for  the  motion  of  compacted  material,  and  hence  of  the  wall,  returns  us  to  the  elastic-plastic 
case  (Sections  2.3  and  2.4)  -  and  again  material  freezes  just  before  it  reverses  direction,  effectively 

stabilizing  the  cavity. 

The  main  effects  of  compaction,  then,  are  (a)  to  inject  a  second  dissipative  process  into  cavity 
formation,  thereby  causing  more  rapid  decay  of  the  shock,  and  (more  important)  (b)  to  enlarge  the 
ratio  of  wavespeed  in  disturbed  material  to  that  of  the  shock.  The  latter  matter  has  great  practical 
import  because  it  lends  physical  sense  to  the  treatment  of  shocked  material  as  rigid-plastic  (i.e., 
incompressible,  but  with  finite  shear  strength),  while  reserving  to  the  undisturbed  medium  ordinary 
compressibility  in  shock  loading -a  locking  solid  with  strength.  This  in  turn  has  practical  impact 
because  the  problem  of  spherical  motion  in  such  a  medium,  when  driven  by  radius-dependent  pressure 
at  the  cavity- wall  (an  arbitrarily  specified  pressure),  has  been  solved  exactly  (Trulio,  1992) -even 
allowing  strength  to  vary  with  mean  stress  (as  it  often  does  in  compactible  media). 
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SECTION  3 

APPROACHES  TO  THE  PROBLEM,  AND  THE  ROUTE  CHOSEN 
3.1  CRATERS  FROM  CAVITIES. 

Evidently,  a  rather  detailed  answer  can  be  given  to  the  question  of  how  cavities  form  around  contained 
explosions,  but  what  has  that  to  do  with  the  craters  that  explosions  produce?  Well -the  two  are 
linked  in  a  general  way  by  the  mechanisms  of  flow  and  compaction,  which  play  a  major  role  in  the 
formation  of  both  cavities  and  craters.  Moreover,  the  linkage  becomes  specific  and  quantitative  in  that 
the  field  around  a  spherical  cavity  automatically  presents  us  with  an  instance  of  ground-surface 
loading  —  one  with  important  similarities  to  the  loads  imposed  by  surface  bursts.  In  a  spherical  field 
around  a  cavity,  the  loading  in  question  is  that  found  in  the  medium  on  any  plane  -  a  horizontal  plane 
9,  say  —  through  the  burst-point;  the  hemispherical  portion  of  the  cavity  wall  below  '5P  completes  the 
surface,  and  the  pressure  exerted  on  it  by  cavity  gases  is  clearly  relevant  to  surface-burst  loading  of  the 
ground.  As  noted  above,  however  (Section  2.1),  the  hoop  stress  imposed  in  plane  on  the  medium 
would  be  tensile  and  hence  qualitatively  (i.e.,  grossly)  different  from  surface-burst  loading,  if  the 
medium  behaved  elastically. 

The  exact  fields  at  our  disposal  are  not  limited  to  elasticity,  though;  they  include  shear  failure  and 
compaction,  which  is  fortunate  even  for  understanding  cavity  growth,  since  models  incorporating  those 
processes  better  represent  geo-materials  than  does  elasticity  alone.  And,  in  plane  <5>,  hoop  stress 
presents  a  quite  different  picture  for  elasto-plastic  media  than  for  elastic  ones;  even  in  quasistatic 
spherical  fields  (as  also  noted  above;  Section  2.3),  hoop  stress  is  compressive  almost  to  the  elastic- 
plastic  boundary  (85%  of  the  way  out).  In  locking  solids  (with  or  without  compaction),  it  remains 
compressive  to  at  least  the  elastic-plastic  boundary  (Section  4  plots  it  in  specific  cases).  Moreover,  the 
material  freezes  (above)  when  the  outgoing  shock  reaches  that  boundary;  i.e.,  since  material  beyond  the 
boundary  resists  shear  elastically,  elements  of  material  thereafter  execute  only  damped  displacement 
oscillations  of  small  amplitude  about  their  final  radii.  Up  to  the  time  of  freezing,  a  loading  wave 
sweeps  outward  over  plane  ?P  imposing  compressive  normal  tractions  thereon,  while  gaseous  products  of 
the  explosion  press  on  the  wall  of  the  expanding  crater  that  these  tractions  produce  (half  the  cavity  of 
the  contained-burst  field  whence  the  tractions  come)  -  all  in  qualitative  agreement  with  surface-burst 
loading. 

The  field  below  ?P  also  differs  strikingly  from  that  of  a  surface  burst,  by  presenting  such  major 
simplifications  as  a  hemispherical  crater  (not  merely  bowl-shaped)  and  a  just-sonic  blast  wave; 
noteworthy  too  is  a  discontinuity  in  loading  pressure  of  amount  Y  [for  simple  materials  (Appendix  B.l) 
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during  flow]  on  moving,  at  the  crater’s  edge,  from  the  crater-wall  (vertical)  to  flat  ground  (horizontal). 
Yet,  the  surface  tractions  are  of  the  right  kind:  compressive,  spreading  and  diminishing  in  amplitude 
both  with  time  and  with  distance  from  the  burst  point.  The  consequent  fractions  of  crater  volume  due 
to  compaction,  flow  and  elastic  volume-transport  therefore  warrant  serious  attention;  they  could  prove 
representative  of  exact  solutions  for  true  surface  bursts,  and  of  the  trends  exhibited  therein. 
Encouraging  of  such  a  prospect  is  the  fact  that  (Trulio;  1984),  for  bursts  at  and  near  the  flat  surfaces  of 
homogeneous  halfspaces,  motion  is  almost  spherical  in  a  cone  of  angle  60°,  at  least,  that  opens 
downward  from  its  vertex  directly  above  shot-point  and  grazes  the  crater’s  rim.  Indeed,  the  present 
line  of  attack  is  less  ambitious  than  attempts  via  “equivalent  yield”  (which  the  approach  here 
resembles)  to  relate  surface  bursts  to  contained  ones.  That  is,  since  the  volume-fractions  at  issue  are 
independent  of  yield,  no  absolute  correspondence  need  be  made  between  the  two  types  of  burst. 

As  to  calculating  the  volume-fractions,  some  large  points  remain  to  be  covered.  The  key:  Given  the 
state  of  the  medium  (out  to  the  elastic-plastic  boundary)  after  it  freezes,  we  assume  that  almost  the 
same  state  would  have  been  reached  had  the  medium’s  compressibility  been  at  the  low  levels 
appropriate  to  elastic  unloading  after  compaction.  Pertinent  compressibilities,  where  needed,  come 
ultimately  from  lab  measurements  on  materials  of  interest.  Independently  of  those,  but  given  that 
deformation  is  purely  elastic  outside  the  elastic-plastic  boundary,  the  volume  transported  elastically  to 
great  distances  equals  that  between  the  final  and  initial  radii  of  material  at  that  boundary  (a  spherical 
shell).  Likewise,  without  reference  to  actual  non-zero  compressibilities,  the  exact  locking-solid  solution 
provides  the  volume  compacted  out  of  material  between  the  cavity-wall  and  the  elastic-plastic 
boundary.  The  third  piece  of  the  crater  -  the  part  added  to  its  volume  by  flow  -  could  be  viewed  as 
the  volume  transported  elastically  plus  that  taken  up  in  compression  by  locked-in  stresses  (above). 
Having  already  accounted  for  the  former,  however,  flow  is  associated  here  only  with  locked-in  stress. 
The  volume  due  to  flow  then  follows  from  the  mean-stress  distribution  out  to  the  elastic-plastic 
boundary,  using  the  medium’s  non-zero  compressibility. 

Finally,  one  takes  account  of  a  patent  difference  between  surface  and  contained  bursts,  by  unloading 
the  ground  surface  to  zero  normal  stress  (.1  MPa,  but  who’s  counting).  This  last  portion  of  the 
calculation,  the  only  one  that  truly  involves  two  spatial  dimensions,  has  two  parts.  First,  the  cavity 
wall  is  relieved  elastically  of  its  uniform  pressure  by  means  of  spherical  deformation  (by  superposition, 
using  the  known  field  that  results  elastically  from  such  relief),  which  alters  the  radial  distribution  of 
normal  (hoop)  stress  over  the  flat  part  of  the  surface.  The  latter  adjusted  stresses  are  then  relieved 
elastically,  also  by  superposition,  but  using  the  known  two-dimensional  field  produced  by  that  relief 
(Love,  1944,  p.  192).  In  the  resultant  field,  the  elastic-plastic  boundary  will  have  moved,  becoming 
aspherical,  and  the  volume  between  its  former  and  new  loci  must  be  added  to  the  prior  estimate  of  the 
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volume  transported  elastically  (or  the  entire  calculation  of  said  volume  can  await  fixing  the  final 
locus).  Similarly,  the  mean-stress  distribution  between  the  cavity  wall  (aspherical  in  the  end)  and 
elastic-plastic  boundary  will  have  changed  too,  requiring  an  addition  to  the  prior  estimate  of  volume- 
due-to-flow  (or  calculation  of  that  volume  from  the  final  field). 

3.2  RELATED  EXACT  FIELDS  IN  ELASTIC  MEDIA  AND  INCOMPRESSIBLE  FLUIDS. 

The  process  outlined  plainly  provides  no  more  than  approximate  final  fields  for  surface  bursts,  and 
hence  approximate  crater-volume  budgets  for  those  axisymmetric  fields  (though  the  volume-fractions 
desired  may  well  be  more  accurate  than  field- values  at  individual  points).  It  does  so,  however,  for 
media  that  deform  inelastically,  both  hydrostatically  and  in  shear,  as  well  as  elastically.  No  exact 
axisymmetric  fields  are  known  for  such  media.  Rather,  the  roster  of  exact  axisymmetric  fields  is 
limited  to  linear  elastic  materials,  using  linearized  equations  of  motion,  and  to  incompressible  fluids. 
Those  had  been  proposed  nonetheless  (for  this  program)  as  possible  sources  of  estimates  of  elastic- 
transport  volumes.  One  could  learn,  for  example,  where  the  volume  of  a  dent  goes  when  a  localized 
load,  impulsively  delivered,  is  impressed  on  a  linear  elastic  halfspace;  as  one  item  in  question, 
permanent  non-zero  volume  strains  might  develop  near  the  loaded  surface  (if  the  load  is  not  removed), 
or  a  part  of  that  surface  may  be  raised  above  its  initial  plane.  Further,  one  would  see  how  compression 
varies  in  the  outgoing  shock  with  direction  from  the  central  point  of  loading.  At  a  more  basic  level, 
one  might  find  that  a  local  load,  even  though  it  delivers  infinite  impulse  (if  maintained  indefinitely), 
will  not  displace  the  entire  halfspace  by  a  finite  distance  [the  exact  solution  for  quasistatic  loading  (loc. 
cit.)  seems  to  imply  that  it  will  not];  if  the  local  load  is  relieved,  so  that  the  impulse  delivered  is  finite, 
then  displacement  of  the  entire  halfspace  does  not  seem  possible  [in  the  simple  case  of  uniform  loading 
(above),  it  takes  infinite  impulse  to  cause  finite  permanent  displacement  of  the  entire  surface,  so  that 
finite  impulse  should  impart  no  permanent  displacement  to  it].  Moreover,  questions  like  these  would 
be  answered  for  the  quintessential  simple  solid  —  but  at  the  price  of  a  small-strain,  small-displacement 
requirement  that  rules  out  indentations  with  crater-like  shape. 

On  a  halfspace  of  incompressible  fluid,  surface-burst  loading  would  produce  a  crater  which,  under 
gravity,  would  eventually  disappear.  Without  gravity,  it  might  stabilize  at  finite  size  — but  if  so, 
would  surely  be  huge  relative  to  actual  craters.  That  point  aside,  making  room  for  the  crater  would 
call,  in  equal  volume,  either  for  a  net  outflow  at  great  distances  (across  a  vast  spherical  surface  about 
the  burst-point,  say),  or  for  the  elevation  of  some  or  all  of  the  medium’s  surface  outside  the  crater; 
changes  in  material  density  are  excluded  by  definition.  The  medium  has  charm  nonetheless,  because 
calculating  exact  flow-fields  for  it,  large  strains  and  displacements  included ,  is  a  tractable  problem  (but 
not  simple;  solving  Laplace’s  equation  is  the  chief  task).  On  the  other  hand,  neither  option  (elastic 
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materials  or  incompressible  fluids)  permits  compaction;  further,  flow  too  is  excluded  from  elastic 
halfspaces,  as  is  elastic  deformation  from  incompressible  fluids -each  being  a  major  omission  from  the 
mechanics  of  crater  formation. 

3.3  DISCRETE  METHODS,  AND  THE  APPROACH  CHOSEN. 

Massive  computation  based  on  discrete  equations  of  motion -a  favored  option  in  the  computer 
age -was  relegated  from  the  outset  to  no  more  than  a  supporting  role  (mainly  to  be  played  later,  for 
purposes  of  comparison).  To  be  sure,  discrete  methods  have  the  advantage  of  readily  admitting  in  a 
material  the  modes  of  behavior  noted  (and  others,  to  the  point  of  diminishing  returns).  The  aim  here, 
however,  is  to  cover  a  wide  range  of  those  properties,  and  each  choice  of  the  set  of  them  gives  rise  to  a 
distinct  cratering  field.  Calculating  one  such  axisymmetric  field  is  expensive,  and  the  result  becomes 
exact  only  in  the  limit  of  infinite  expense.  In  practice,  that  is,  discrete  methods  always  give  inexact 
answers  to  the  specific  problems  of  motion  posed,  and  their  cost  is  compounded  when  one  seeks  high 
accuracy  at  late  times  and/or  large  distances. 

In  particular,  for  axisymmetric  fields,  the  number  of  arithmetic  operations  must  grow  as  the  cube  of 
final  problem  time  if  we  are  to  have  much  hope  of  controlling  discretization  error.  Even  then,  error- 
control  is  less  than  certain:  It  has  not  been  shown  that  discretization  error  can  be  kept  in  check  as  a 
cratering  problem  is  run  on  and  on,  and  to  make  such  an  evaluation  with  confidence  would  be 
expensive  and  time-consuming.  Yet,  for  present  purposes,  a  numerical  field  without  a  demonstrable 
bound  (at  least)  on  its  error  would  be  worthless.  Likewise,  so  are  such  tricks  as  computing  the  near 
field  using  discrete  cells  of  the  usual  size,  and  then  equating  the  volume  transported  elastically  to  the 
final  computed  crater  volume  minus  the  sum  of  “permanent”  volume  changes  over  all  cells  of  the 
discrete  mesh.  For  one  thing,  without  artificial  damping  (explicit  or  implicit  in  the  discrete  equations 
of  motion),  the  mesh-nodes  for  such  a  problem  oscillate  indefinitely,  and  it  is  not  known  how  much  the 
final  axisymmetric  density  field,  damped,  differs  from  the  undamped  one  sought.  The  displacement  of 
material  at  a  discrete  elastic-plastic  boundary  is  subject  as  well  to  that  uncertainty. 

In  sum:  (i)  Cost,  of  itself,  sufficed  to  tip  the  scales  against  massive  computation,  (ii)  Either  course 
(heavy  computing,  or  using  exact  solutions  to  related  problems)  yields  only  approximate  results,  but, 
with  regard  to  understanding  the  defects  of  the  crater-volume  budgets  obtained,  clearcut  physical 
compromises  are  much  to  be  preferred  to  the  obscure  (and,  mathematically,  no  less  complex)  effects  of 
discretization  error.  Thus,  early  in  the  program,  a  choice  had  to  be  made  between  exact  solutions  that 
take  account  of  just  one  of  the  three  physical  processes  most  basic  to  cratering,  and  related  exact 
solutions  which  include  all  three.  Given  the  objective  of  determining  what  these  very  processes  might 
contribute  to  final  crater  volume,  and  almost  no  prior  estimates  thereof,  we  chose  to  employ  model- 
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materials  capable  of  flow,  compaction  and  elastic  deformation  -  and  compromised  on  completeness  and 
exactness  in  other  aspects  of  the  problem  (chiefly  ground  loading  and  the  passage  from  locked  solids  to 
stiff  elast  ic  ones).  In  short,  the  approach  outlined  in  Section  3.1  was  adopted. 
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SECTION  4 

THE  FIELDS  CALCULATED:  MEDIA,  PROCEDURE,  AND  RESULTS 


4.1  CONSTITUTIVE  EQUATIONS. 

In  calculations,  including  those  reported  here,  models  define  materials.  The  models  take  the  form  of 
equations  (“constitutive  equations”)  which,  ultimately,  relate  stress  and  strain  increments  in  any  given 
element  of  material.  Each  such  equation-set  defines  a  class  of  substances.  Individual  members  of  the 
class  are  specified  by  the  values  assigned  to  all  the  arbitrary  (within  limits)  constants  contained  in  the 
equation-set  (material  properties  “parameters”).  Comprising  the  constitutive  equations  used  m  this 
work  are:  (i)  a  curve  of  stress  vs.  strain  (uniaxial)  for  shock  loading,  from  whose  states  all  subsequent 
stress-strain  changes  take  place  at  constant  volume  (reloading  to  higher  mean  stresses  than  m  the 
initial  shock,  if  that  occurs,  is  then  easy  to  handle);  (ii)  the  requirement  that,  in  shock-states,  for  which 
<rr  and  ae  are  denoted  <r*  and  o%  ,  \<rr-<re\  equal  the  smaller  of  Y  and  ar{ and  that 
otherwise  (iii)  |<rr  -  ag\  equal  the  smaller  of  Y  and  |<  -o%-  6Ge„|,  with  <rr  -  tr0=*Br  -<rse-  6Gee  when 
Y  is  the  larger  number.  In  calculating  elastic  transport  of  volume  outward  from  the  elastic-plastic 
boundary,  material  beyond  the  boundary  obeys  the  equations  of  isotropic  linear  elasticity,  which  (in 
this  context)  reduce  to  <rr=A  A  +2Ger  and  ae=  A  A  +2Ge„,  where  A  =  er+2e0;  the  same  equations  are 
used  to  compute  the  elastic  compression  that  locked-in  stresses  would  cause.  The  model’s  parameters 
are  then  the  constants  Y,  v,  G,  A,  and  those  of  the  loading  curve,  for  which  two  different  algebraic 
forms  were  used:  <r*/«T0=(e */e0)a  [“power-law”]  and  </<r0=a(e*/e0)  +b«/e0)2+b«/e0)3  [“cubic”], 
whose  parameters  <r0,  e0,  a,  a,  b  and  c  occur  in  just  four  independent  combinations  (one  in  the  power- 
law  equation  for  <r‘/<r0  and  three  in  the  cubic).  Compared  to  most  geo-material-models  now  in  use, 
this  one  is  quite  simple,  as  befits  an  exploratory  program  of  the  present  kind.  An  equation  for  pressure 
(radial  stress)  at  the  wall  of  the  cavity  completes  these  problems’  specification;  said  pressure  was  made 
a  function  of  cavity  volume  (hence  radius)  by  means  of  accepted  equations  of  state  for  detonation 
products  (homogenized;  Trulio  1992,  Section  3.2). 


The  burden  and  results  of  the  PEACE  program  having  prompted  this  work,  the  parameters  noted  were 
assigned  values  derived  mainly  from  lab  tests  of  either  Pacific  coral  or  (lacking  data  for  coral) 
limestone.  Thus,  in  the  four  principal  media  for  which  volume  budgets  have  been  obtained  under  the 
subject  contract,  Y  was  set  at  10  MPa -twice  or  thrice  that  actually  measured  for  Pacific  coral 
(Borschel,  Klauber  and  Early,  1986,  pp  81  k  93)  because  the  present  application  spans  several  times 
the  depth  whence  came  the  samples  tested,  and,  deeper  down,  the  medium  may  be  stronger.  Further, 
v  was  assigned  the  value  .21  (Borschel,  Klauber  and  Early,  1986,  p  129)  in  all  four;  in  one  of  them, 
however,  v  was  also  set  to  0  (A=0)  and  -1  (the  lowest  value  possible  elastically)  with  all  other 
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parameters  held  fixed.  Three  of  the  media  had  power-law  loading  curves  (with  exponents  2,  8.4  and  1); 
the  fourth  (the  one  with  which  iz-values  of  .21,  0  and  -1  were  used)  had  a  cubic  curve  of  loading.  At 
low  stress  levels,  the  latter  curve  fits  uniaxial  stress-strain  data  for  Pacific  coral  (Blouin  and  Timian, 
1986,  p  124);  it  also  fits  selected  stress-strain  points  well  above  the  limit  of  those  data  (.17  GPa  of 
axial  stress),  conforming  roughly  to  Hugoniot  measurements  for  Indiana  limestone  at  10  GPa  (Murri, 
Grady  and  Mahrer,  1975).  The  loading  curve  so  defined  represents  fairly  well  the  data  cited  over  the 
range  of  shock  stress  met  with  here:  0  <  a\  <  14.5  GPa.  The  power-law  curve  with  a=2  is  a  fit  to  the 
same  lab-test  data  for  axial  stresses  up  to  ~  .17  G P a  (though  used  to  much  higher  stresses).  The 
curve  with  a— 1  passes  through  the  10-GPa  point,  but  is  intended  mainly  to  set  a  plausible  minimum 
on  energy  dissipation  in  shock  loading  (zero)  because  the  Rayleigh  line  and  loading  curve  then  coincide; 
a=8.4  comes  close  to  maximizing  such  dissipation,  besides  which,  with  o=8.4,  a  close  fit  is  obtained  to 
the  measurements  of  axial  stress  at  10  GPa  (above),  and  to  the  highest-stress  point  (  ~  .17  GPa)  of  one 
test  in  the  set  whose  overall  data  suggested  o=2.  The  constant  <r0  is  equal  to  the  mean  pressure  of 
cavity  gases  when  detonation  ends  (given  a  chemical-explosive  sphere,  tamped  in  the  medium  and 
center-detonated).  Then,  having  specified  o  for  a  power-law  loading  curve,  e0  is  obtained  by  forcing 
the  curve  to  pass  through  a  specific  stress-strain  point  (at  10  GPa  for  0=1;  .17  GPa  for  a  8.4  and  2); 
for  the  cubic,  whose  coefficients  depend  on  the  choice  of  e0,  the  value  chosen  for  it  forces  their  sum  to 
equal  1  when  <r*=<r0,  and  the  curve  must  go  through  (<r0,c0).  Further,  exclusive  use  of  the  explosive 
Com  position  C4  (as  defined  by  a  JWL  equation  of  state;  Lee  and  Finger,  1974)  to  drive  the  medium, 
led  to  the  value  14.5  GPa  for  a0.  In  some  respects,  especially  in  the  shape  of  the  curve  of  uniaxial 
loading  and  the  “Poisson’s  ratio”  v,  a  broad  range  of  materials  was  covered  by  this  pilot  study;  for 
others  (Y  in  particular)  its  coverage  was  highly  limited  or  minimal. 

4.2  THE  SPECIFIC  PROBLEMS  FOR  WHICH  EXACT  SOLUTIONS  WERE  FOUND. 

Table  4-1  below  lists  the  combinations  of  parameter  and  loading-curve  choices  defining  the  media  for 
which  we  computed  crater-volume  budgets.  As  explained  above  (Section  3),  exact  fields  driven 
explosively  into  locking  solids  form  the  basis  for  the  various  budgets,  and  all  the  properties  defining 
those  solids  appear  in  the  table  (each  row  thereof,  and  hence  a  solid,  is  labeled  in  the  table  s  first 
column  simply  by  the  number  of  the  locking-solid  problem  solved  in  evaluating,  for  that  solid,  the 
exact  solution).  The  “Shear  Modulus,  G,”  listed  in  the  last  column,  does  not  actually  enter  the  exact 
locking-solid  solution,  however.  True,  once  shocked,  a  thin  spherical  shell  of  material  continues  to  fail 
in  shear  if  and  only  if,  on  assuming  its  subsequent  changes  in  strain  to  be  elastic,  one  obtains  principal- 
stress  differences  >  Y.  Since  all  deformation  takes  place  at  constant  volume-strain  A  (Section  4.1), 
this  means  that  shear  failure  (trr  —  <r^=Y;  rigid  plasticity)  continues  if  and  only  if  \<rr  —  —  6Gfg|  >  Y 

as  the  shock  runs  past  the  shell.  Now,  is  zero  in  just-shocked  material  because  displacement  cannot 
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Table  4-1.  Parameters  used  in  material  models. 


Problem 

Number 

Loading 

Curve 

€o 

a 

Poisson’s 
Ratio,  v 

9 

Power- law 

1.68939556 

2.0 

0.  21 

10 

Power-law 

0.31368232 

8.4 

0.  21 

11 

Power-law 

0.36361471 

1.0 

0.  21 

12 

Cubic 

0.29678248 

* 

0.  21 

12a 

Cubic 

0.29678248 

* 

0.  00 

12b 

Cubic 

0.29678248 

* 

-1.  00 

Composition  C4 
14.5445838  GPa 
10.0  MPa 
2.0  GPa 

*  Coefficents  of  cubic 
a  =  .1081465913 
b  =  .5664047702 
c  =  .3254486385 


For  all  problems: 
Explosive: 

Y: 

G: 
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jump  in  an  instant.  But,  because  the  shock  imparts  outward  velocity  to  the  shell,  e0  becomes 
increasingly  tensile  (<0)  thereafter  so  that,  with  o*T-o%  =Y,  we  find  6Ge^>0 

_  (j8q  —  6G€^=<Tr  —  <tq  for  rigid-e/asfoc  deformation).  Hence,  the  shell  continues  to  fail  in  shear.  By 
the  same  token,  if  cr*  -  were  to  fall  below  Y  beyond  some  radius  (due,  say,  to  decay  of  shock-stress 
levels  as  the  shock  moves  outward),  ar  —  (Tq  would  return  at  once  to  the  value  Y.  Thus,  by  driving 
material  outward,  the  shock  also  forces  it  to  fail  in  shear;  moreover,  outward-moving  material  keeps  on 
failing  in  shear.  In  turn,  since  compressive  shock  stress  causes  outward  motion,  the  medium  moves  in 
shear-failure  states  until  the  shock’s  strength  decays  to  zero -which,  given  shear  failure,  happens  at  a 
finite  radius  [Trulio,  loc.  cit.;  cohesion  (in  this  case  Y)  enters  the  shock  stress  equation  as  a  direct 
reduction  of  cavity  pressure  because  flow  at  fixed  shear  stress  Y/2  need  only  occur  in  a  finite  mass  in 
order  to  dissipate  the  explosive’s  finite  energy]. 

To  introduce  elastic  volume-transport  into  such  an  exact  solution,  the  medium  is  treated  as  linear- 
elastic  beyond  an  elastic- plastic  boundary  defined  as  follows:  For  each  position  of  the  shock  (out  to 
that  boundary)  the  radial  velocity  of  a  just-shocked  shell  is  known,  and  cannot  jump  discontinuously  at 
any  radius  (because  it  would  then  have  to  jump  throughout  the  mass  of  locked  solid,  which  would  take 
infinite  force).  If  a  shock  were  driven  by  locked  solid  into  elastic  material  just  beyond  a  given  radius, 
the  radial  velocity  found  in  the  locked  material  just  shocked  would  therefore  be  imparted  to  the 
adjacent  elastic  shell.  Values  having  been  assigned  to  that  shell’s  elastic  moduli  and  density,  its 
longitudinal  wavespeed  is  also  known  —  and  with  it,  radial  stress  in  the  shell’s  shock-state  (the  product 
of  density,  wavespeed  and  radial  material  velocity),  hoop  stress,  and  their  difference  ar  —  (Tq  [radial 
stress  times  (1-2i/)/(1-i/),  since  the  shock  state  is  one  of  uniaxial  strain].  The  elastic-plastic  boundary 
lies  (before  displacement  develops  there)  at  the  radius  for  which  —  so  computed  (i.e.,  in  the 
elastic  shell),  equals  Y. 

Since  the  moduli  define  a  material  much  stiffer  in  compression  than  the  loading  curve  requires  (at  low 
stresses),  radial  stress  jumps  to  higher  levels  than  in  the  shock  (through  locking  solid)  incident  upon 
the  boundary.  That  compressive  jump  is  felt  at  once  throughout  locked  solid  as  a  similar  radius- 
dependent  jump  (zero  only  at  the  cavity  wall);  velocity  must  be  continuous,  but  not  acceleration. 
Owing  to  the  great  stiffness  of  the  elastic  medium,  particularly  in  shear,  the  system  now  freezes 
(Section  2.3).  That  is,  small  tensile  changes  in  hoop  strain  (hence  in  radial  displacement)  cause 
relatively  large  tensile  changes  in  hoop  stress  in  the  outgoing  elastic  wave,  soon  bringing  outward 
motion  of  the  elastic-plastic  boundary  to  a  halt,  and  the  boundary  starts  to  move  inward.  Now, 
however,  hoop  strain  becomes  more  compressive  in  locked  solid,  causing  <rr  —  <t$  to  fall  below  Y,  so  that 
(finally)  rigid  -elastic  deformation  of  locked  material  occurs;  the  large  shear  modulus  for  such 
deformation  restricts  the  oscillation  to  virtually  negligible  displacements.  The  exact  solution  extends  to 
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this  spherical  system  as  well:  locking  solid,  elastic  (in  shear)  and  plastic,  surrounded  by  an  elastic 
medium.  Volume  transported  elastically  issues  from  it;  a  close  approximation  thereto  comes  from 
applying  the  shock  stress  incident  on  the  elastic-plastic  boundary  to  the  outside  elastic  medium,  as  a 
quasistatic  load. 

In  the  shear-failure  region,  where  larger  compressive  stresses  are  found  than  outside  its  elastic-plastic 
boundary,  elastic  moduli  for  unloading/reloading  to  vary  with  shock  stress.  Specifically,  data  from  UX 
tests  of  Pacific  coral  show  nearly  elastic  unloading/reloading  from  a  given  loading  state,  but  with 
moduli  M  and  G  that  depend  appreciably  on  the  loading  state  chosen  (Borschel,  Klauber  and  Early, 
loc.  cit.).  Nonetheless,  G  and  the  Poisson’s  ratio  vu  for  unloading/reloading  were  assigned  constant 
values  (2  GPa  and  .21,  respectively,  found  at  the  higher  stresses  of  UX  measurement);  from  the  usual 
definitions  of  elastic  moduli,  G  and  vu  imply  a  constrained  modulus  M,  and  also  a  bulk  modulus  K 
(=M-iG).  For  a  given  spherical  shell  of  material,  the  locked-in  elastic  volume-strain  (the  change  in  its 
volume-strain  that  would  occur  on  unloading  to  overburden)  equals  its  locked-in  mean  stress  <r 
[=i(<Tr+2<rfl)],  divided  by  K;  thus,  the  product  of  (f  /K  with  the  shell’s  final  volume  gives  its 
contribution  to  locked-in  volume.  Summing  over  all  shells  in  the  shear  failure  region  [i.e.,  evaluating 
/(<?  /K)4flT2  dr  across  the  shear-failure  region],  we  obtain  the  total  volume  locked  into  that  region  — the 
contribution  of  flow  to  the  crater’s  volume. 

Formulas  for  the  change  in  the  field  due  to  unloading  of  the  ground  surface  are  those  of  strict  linear 
isotropic  elasticity,  and  hence  contain  just  two  constant  moduli  (bulk  and  shear  moduli,  say).  Long 
since  published  (Love,  loc.  cit.),  they  yield  the  following  expressions  for  the  vertical  (uj  and 
horizontal-radial  (uR)  components  of  displacement  at  an  arbitrary  vector-position  r  in  the  halfspace: 

r  *  dx^d  * 

4ttGu  =  2(1-1/)^  ~zd<j>/dz  ;  <£=  ^„(x',y')  -nMr  C4-1) 

J  J  GS  11  1 1 


4ttGur  =  —  z5<^/5R  — (1-2j/)5x/9R 


dx/dR  = 


n(x',y ')  [R-R'Cos(y>-y?<)]  d  , 
3  |r-r'|(z+|r-r'|) 


where  integration  extends  over  the  part  of  the  ground-surface  plane  [the  (x,y)-plane  of  a  Cartesian 
(x,y,z)  frame]  loaded  with  non-zero  normal  stress  <r„(x,y)  (hoop  stress  in  the  exact  spherical  solution, 
after  unloading  the  cavity  wall;  Section  2.1);  also,  z  denotes  vertical  position,  increasing  with  depth, 
x',y',0  are  the  position-components  of  an  arbitrary  point  on  the  ground-plane.  Points  r  and  r'  have 


cylindrical  coordinates  R,<^  and  R' (x— R  Cos^>,  y— R  Siny>;  x  —  R  Cosy* ,  y  R  Sin<£> ). 


Here  «rn(x',y'), 


a  function  only  of  R',  is  the  normal  stress  at  r'  when  the  medium  freezes,  so  that  the  region  GS 


becomes  the  circle  (about  shot-point)  whose  radius  equals  the  elastic-plastic  radius  r£p;  passage  to 
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cylindrical  coordinates  is  therefore  convenient,  whereupon  R'dR'dy?'  replaces  dx'dy ,  0  <  R  <  rep  and 
0  <  y>'  <  2ir,  and  <j>  and  8X/d R  (hence  uz  and  uR)  are  independent  of  <p  [also, 

|r-r,|2=z2+R2+R'2-2RR'Cos(y>-v?')]- 
4.3  RESULTS. 

At  about  the  time  it  freezes  in  our  exact  spherical  field,  medium  12  (a  locking  solid  with  strength, 
Table  4-1)  exhibits  the  hoop-stress  distribution  plotted  as  a  function  of  radius  in  Figure  4-1.  That 
stress  (as  it  evolves  with  time),  if  imposed  normally  on  a  flat  halfspace  of  the  same  material  (except  for 
a  hemispherical  indentation  loaded  with  the  cavity’s  pressure),  produces  every  detail  of  the  field  found 
spherically.  As  a  ground-surface  load,  it  is  therefore  compared  in  the  figure  with  two  others  that  come, 
respectively,  from  computer  simulations  of  conventional  and  nuclear  surface  bursts  (zero  height-of- 
burst)  —  the  Middle  Gust  4  event  (Ito,  1995)  and  an  unspecified  one  (Brode,  1964)  Since  no  known 
scaling  law  relates  the  three  fields,  the  general  idea  of  crater-volume  scaling  (Cooper,  1976)  was 
employed  in  selecting  the  overpressure  profiles  shown.  Specifically,  the  profiles  are  found  in  each  of  the 
two  latter  fields  at  about  the  time  when  the  loading  wave  has  moved  outward  by  the  number  of  cavity 
radii  seen  in  the  hoop-stress  profile  [the  cube  of  cavity  radius  at  a  given  time  being  defined  as  then- 
current  cavity- volume  x  (.75 / tt)]-  Further,  differences  among  the  three  media  were  taken  into  account 
at  the  outset  by  applying  to  current  crater  volume  a  crude  empirical  factor  that  allows  crater  size  to 
vary  with  medium  (soil,  soft  rock  and  hard  rock;  Cooper,  1976). 

The  shapes  of  the  surface-burst  profiles  are  plainly  much  more  alike  than  either  is  to  the  hoop-stress 
profile,  displaying  peak  values  not  near  the  crater  but  farthest  away  from  it  (i.e.,  at  the  front  of  the 
loading  wave).  That  fact  evidences  the  different  physical  objects  -  airblast  front  vs.  locked-in 
stress  -  presented  by  the  two  types  of  peak.  Peak  loading  stresses  are  about  7  times  as  large  in 
medium  12  as  for  the  surface  bursts,  and,  in  the  more  important  matter  of  vertical  force  exerted  on  the 
ground,  the  ratio  is  ~  10  —  but,  at  the  wavefront,  the  stress-ratio  is  <j.  Thus,  the  two  types  of  ground- 
loading  might  give  rise  to  similar  crater-volume  contributions,  percentagewise,  from  elastic  transport. 
The  higher  stresses  in  medium  12,  however,  argue  that  the  present  fields  overstate  the  contribution  to 
crater  volume  from  compression,  elastic  or  not.  To  answer  such  questions  requires  a  better  match  to 
surface-burst  loading  than  is  provided  by  the  hoop-stress  profiles  from  contained  bursts.  Such 
improvement  can  result  from  varying  the  cavity-radius/wall-pressure  relation  by  which  our  exact  fields 
are  driven.  The  exact  solution  offers  great  leeway  in  setting  that  relation,  which  was  in  fact  varied  for 
medium  12  to  produce  the  more  satisfactory  match  to  surface-burst  loading  seen  at  the  bottom  of 
Figure  4-1;  to  obtain  the  hoop-stress  curve  shown  there  required  dictating  far  more  rapid  decay  of 
cavity-wall  pressure  than  that  due  to  cavity  gases  in  a  contained  burst -aping,  crudely,  the  expansion 
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Figure  4-1.  Comparison  of  normal  stress  loads  on  the  ground  surface  from:  This  paper  [IIoop  Stress 

(Medium  12)],  a  nuclear  surface  burst  (Brode)  and  a  conventional  burst  (TRT).  For  graph 
I,  the  pressure  of  gaseous  detonation  products  drives  the  field  whose  Hoop  Stress  is  plotted. 
Graph  II  employs  a  faster-falling  relation  between  cavity  pressure  and  radius,  contrived  to 
make  Hoop  Stress  match  more  closely  the  TRT  and  Brode  ground-surface  loads. 
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of  gaseous  detonation  products  into  the  atmosphere  during  surface  bursts. 

As  for  the  three  contributions  themselves  (locked-in  elastic  volume,  elastic  transport  and  compaction), 
Table  4-2  lists  the  fractions  of  final  crater  volume  due  to  each,  for  the  six  media  of  Table  4-1.  The 
fraction  due  to  elastic  transport  ranges  from  about  |  to  ±  a  fairly  narrow  range  considering  the  large 
variations  in  shock-loading-curve  and  Poisson’s  ratio  covered  by  these  media  (but  Y,  probably  at  least 
as  important  a  material  property  in  determining  the  volume-fractions  of  Table  4-2,  was  not  varied). 
Most  sensitive  to  changes  in  the  medium’s  properties  is  the  portion  of  crater-volume  due  to  its  elastic 
compression  by  locked-in  stresses;  this  volume-fraction  runs  from  barely  significant  (~i)  t0 
substantial  (  ~  |).  Further,  since  the  boundary  condition  at  the  crater-wall  is  one  of  specified  pressure 
P,  and  the  parameter  P/Y  dominates  quasistatic  cavity  formation,  the  volume-fractions  at  issue  might 
span  a  much  wider  range  if  Y  were  varied.  As  matters  stand,  the  compaction-fraction  is  largest  in  all 
cases  (running  from  .44  to  .80).  Given  the  program’s  variety  of  shock-loading  curves  (of  power-law 
type  especially),  that  result  was  not  foreseen;  e.g.,  with  <r*/<70=.  1,  the  fraction  <*  by  which  volume 
irreversibly  decreases  via  shocking  is  given  by  e*/e0=.  1  when  a=l,  and  c*/£0=.76  when  <*=8.4.  A 
possible  explanation:  Each  of  our  loading  curves  admits  so  much  compaction  that  even  for  the  least 
compatible  medium  of  the  lot  (a=l),  compaction  must  account  for  a  large  part  of  the  space  made 
available  for  cavity  expansion.  High-porosity  Pacific  coral,  after  all,  was  the  medium  whose  measured 
properties  most  influenced  the  sets  of  parameters  for  which  our  exact  field  was  evaluated.  Indeed,  in 
retrospect,  one  might  have  wished  for  less  emphasis  on  the  Pacific  and  more  on  diversity  of  media. 
Prospectively,  more  relevant  halfspace  loadings  such  as  that  at  the  bottom  of  Figure  4-1  [even  if  found 
in  spherical  fields  by  discrete  methods  (fine  zoning)],  should  be  used  to  assess  the  accuracy  of  the 
accounts  of  cratering  that  massive  computation  yields. 

Notwithstanding  the  restriction  to  just  one  Y-value,  insensitivity  of  the  crater-volume  budget  in 
medium  12  to  a  drastic  change  in  pressure  at  the  cavity  wall  (part  II  of  Table  4-2),  and  to  the 
simulated  ground-surface  loading  that  follows  from  wall-pressure  (Figure  4-1),  must  be  counted  a  major 
result  of  this  study.  In  particular,  the  fraction  of  crater  volume  that  resulted  from  elastic 
transport  —  the  study’s  main  concern  —  changed  from  .33  to  .29.  More  strongly  affected  was  the 
volume-fraction  due  to  residual  compression,  which  dropped  from  .22  to  .14  — a  direction  of  change 
consistent  with  the  much-enhanced  decay  of  wall-pressure  that  caused  it.  The  two  reductions  noted 
entail  an  increase  from  .45  to  .57  in  the  compaction-fraction,  because  the  contributions  made  to  crater 
volume  by  elastic  relief  of  pressure  on  the  hemispherical  crater,  and  on  the  ground-surface  plane, 
proved  negligible.  Relief  of  cavity-wall-pressure,  for  example,  made  a  change  of  only  .15%  in  the 
volume  swept  out  by  the  elastic-plastic  boundary  (i.e.,  in  the  volume  transported  elastically  to  remote 
regions).  This  was  so  despite  the  large  pressure,  relative  to  that  in  a  surface-burst  crater,  remaining  at 
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Table  4-2.  Fractions  of  final  crater  volume  contributed  by  elastic  transport,  compaction  and  residual 
elastic  compression  arising  from  flow.  The  fields  of  part  I  were  driven  by  the  pressure  of 
gaseous  detonation  products.  Part  II  employs  a  faster-falling  relation  between  cavity 
pressure  and  radius,  contrived  to  make  Hoop  Stress  match  more  closely  the  TRT  and 
Brode  ground-surface  loads. 


Medium 

Residual 

Compression 

Elastic 

Transport 

Compaction 

9 

.15 

.31 

.54 

10 

.04 

.16 

.80 

ii 

.26 

.30 

.44 

12 

.22 

.33 

.45 

12a 

.20 

.33 

.47 

12b 

.17 

.32 

.51 

12 

.14 

.29 

.57 
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the  cavity  wall  when  the  outgoing  shock  reached  radius  rep-a  shortcoming  of  the  present  cratering 
model  [to  keep  hoop-stress  positive  meant  that  radial  stress  (wall  pressure)  could  be  no  smaller  than  Y 
(10  MPa  here)].  Yet,  the  negligible  effects  of  surface  stress-relief  also  confirm  a  cornerstone  of  the 
model:  When  shear  failure  stops,  so  does  crater  growth;  the  medium  around  the  crater  does  indeed 

“freeze”  at  that  time. 
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APPENDIX  A 

LINEAR  ELASTICITY:  IMPULSE  AND  SPHERICAL  CAVITY  GROWTH 

In  the  usual  limit  of  small  strain  and  displacement,  motion  is  governed  in  linear  elastic  isotropic  media 
by  two  wave  equations,  one  for  the  irrotational  part  of  material  displacement  (or  velocity),  and  the 
other  for  its  solenoidal  part  (Landau  and  Lifshitz,  1970).  The  former’s  waves  move  with  longitudinal 
speed,  while  those  of  the  latter  equation  move  with  transverse  speed.  In  homogeneous  isotropic 
wholespaces,  however,  spherical  fields  of  motion  are  produced  when  pressure  is  applied  uniformly  to  the 
walls  of  spherical  cavities.  Material  displacement  (or  velocity)  is  then  radially  directed  and,  like  all 
properties  of  spherical  fields,  varies  only  with  radius  (measured,  in  this  case,  from  the  center  of  the 
cavity).  Such  displacement  (or  velocity)  fields  are  curl-free  -  i.e.,  irrotational  (Courant  and  Freidrichs, 
1948).  As  a  result,  material  displacement  (or  velocity)  equals  its  irrotational  part,  and  so  satisfies  a 
wave  equation  with  longitudinal  wavespeed.  Moreover,  since  material  moves  on  radial  lines,  and 
displacement  (or  velocity)  varies  only  with  radius  and  time,  we  need  only  express  the  displacement 
vector  and  the  LaPlacian  V2  in  spherical  coordinates  (Stratton,  1941)  to  find  from  the  wave  equation 
that 

c;232ur/dt2=  r-2(d/c>r)(r2dur/dr)  -  2r~2ur  (A.l) 

where  t,  r,  ur  and  cp  denote  time,  radius,  radial  material  displacement  (or  velocity)  and  longitudinal 
wavespeed,  respectively.  Any  given  plane  9  through  the  center  of  the  cavity  is  orthogonal  to  all 
spheres  about  that  center,  and  can  therefore  be  taken  as  the  surface  on  which  the  spherical  coordinate  9 
is  equal  to  tt/2.  The  stress  component  ae  normal  to  9  is  a  principal  stress  (hoop  stress)  -  and,  by 
definition,  compressive  stress  is  taken  as  positive  herein.  If  the  cavity  wall,  whose  radius  is  a,  is 
subjected  impulsively  to  pressure  P  at  time  t=0  -  a  pressure  held  fixed  thereafter  -  then  Equation  (A.l) 
implies  the  following  on  (Cooper  with  7=0): 

<re  =  <t0(£,t)=(-1P/£3){1  -  [( 1  +2/?£2)Cos(kC 6)  +  (l-2£-2/?£2)|Sin(K<$)]  exp(-f<)}  (A-2) 

Equation  (A.2)  holds  for  t  >  0  and  0  <  r-a  <  cpt,  while  <re=0  if  r-a>cpt.  For  r  <  a  and  t  >  0  we 
obtain  <re= P  (on  9)  by  considering  the  cavity  to  be  full  of  gas  at  pressure  P  [the  force  over  the  circle 
subtended  on  9  by  the  cavity  is  thus  normal  to  9  with  magnitude  ira2P;  more  rigorously,  one  obtains 
the  same  result  for  the  force  on  the  hemispherical  cavity  surface  below  9,  by  summing  (over  that 
surface)  the  tractions  Puw  at  all  its  points  (nw  =  unit  outer  normal  at  a  given  point)].  Also,  the  state 
of  zero  stress  (for  P  as  well  as  in  the  medium)  is  that  of  the  medium  infinitely  far  from  the  cavity 
(assumed  hydrostatic).  Further,  the  variable  C  denotes  r-£+l,  where  the  nondimensional  time  r  and 
radius  4  equal,  respectively,  cpt/a  and  r/a;  for  the  solid  with  shear  modulus  G  and  constrained 
modulus  M,  the  constants  /?,  k  and  6  equal  1-2G/M,  M/G  and  (-1+M/G)1/2,  respectively.  We  note 
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too  that  M =/>c2  (where  p  is  the  material’s  density),  and  that  2G/M— (l-2*/)/(l -v)  if  v  denotes 


Poisson’s  ratio. 


The  force  acting  on  ®P  is  normally  directed  (vertically  downward,  say);  the  amplitude  of  the 
force  is  given  by 

F(t)=  |a(1+r)^2jrrdr  =  27ra2|1+To-<)^  =  ira2P-^a2p||  d£/£2-j^  V’(^CW(-«C)d^2| 

=  ^  +  xa2P  J 1  +  exp (-<)  d£/£2  (A.3) 


where 


HU)  =  (l+2/?e2)Cos««)  +  (l-2^2/?e2)|Sin«5) 


Up  to  any  time  t>0,  the  impulse  delivered  to  the  solid  below  plane  ^P,  also  normally  directed,  has 
amplitude  I('r)  equal  to  J^F(t/)dt/  over  the  time-interval  0  <  t  <  t.  That  is,  with  C  =  ^  — £+1. 

\T{t±+  [1+T'!^exp(_<)df|dr'  =  ln(l+r)+r|1+r^«cp(-«0<if^ 
7ra2(a/cp)  J0l1+r  •>!  «  J  JoJi  ^ 


In  the  (T',£)-plane,  the  same  triangular  region  is  defined  by  {0  <  r'  <  r;  1  <  £  <  l+r'}  as  by 
{1  <  £  <  1  +  t;  £-1  <  r'  <  r},  or  {1  <  £  <  1+r;  0  <  <'  <  r-{+l  =  <}.  It  follows  that 

yfr  =  ln(l+r)+  [ ^  f  ’ "  «p(-0  dr'd£ 

7raz(a/c  )  J  t  J  £-i  £ 

f  ..  r  C 

=  ln(l+r)+  X  V’(^C')  exp(-K<')  d£'d£  (A.6) 

J  l  £  J  o 

This  last  integral  with  respect  to  £'  is  elementary;  noting  that  62=-1+2/k,  the  indefinite  integrals  of 
exp(-«C)  Cos«'$)  and  exp(-<')  Sin«'tf)  equal  the  functions  -iexp(-K<0[Cos(«C'<5H  Sin(«C^)] 
and  _lexp(-KC')[5  Cos(«:C^)+Sin(«;C,6)],  respectively.  Hence: 


I(r)/P 

7ra2(a/c 


>  f1+r  df 

p=,n(1+r)+L  W~< 

< 


(l+2/?£2)+(l-2£-2/?£2) 

-  (l+2/?£2)  exp(-KC)[Cos(«C^)-^  Sin(«£6)] 

-  (l-2£-2/?£2)  exp(-«C)[^  Cos(K<6)+Sin(«;<6)] 


f1+T  d£  f  (l+2/?£2)[Cos««M  Sin(ieC«)]  1 
J  j  2£2  |+(l-2£-2/?£2)[6  Cos(K<6)+Sin(/c<<5)]J 


;-!>{ 


£[(l-$)  Cos(«C5)  -  (1+tf)  Sin«$)] 

[6  Cos(K<6)+Sin(fcC<5)]/^+i[(l+<5)  Cos«$)+(l-6)  Sin(«CS)]/£5 


sp(-«C) 
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=  _X_  +  1  /?{(1-6)[Cos(k£6)-6  Sm««)]-0+*)P  Cos(/<6)+Sin«6)]}  exp(-/cC)  I  J  =  0 
+  [r  d<  {[<5  CosK6)+Sin(K<6)]/£-i[(l+6)  Cos(k<6)+(1-6)  Sin«6)]/42}  exp(-<) 

j  0 

=  -l—  +  ±  /3  [(1-2 6-S2)  Cos(k<6)-(1-26+62)  Sin(K^)]  exp(-/c<)  I  J  =  0 
(1+r)  2  s 

[(1+6)  Cos(/c<6)+(l-6)  Sin(/c<6)][  exp(-K<)]/£  |  £  _0 

+  f  dC  {[6  +i«(l+6)2]  Cos«6)  +  [1  +|k(1-62)]  Sin(*;C6)}[exp(-KC)]/£  (A.7) 

Jo 

As  t- *oo,  the  first  line  on  the  right  side  of  Equation  (A.7)  becomes  equal  to  1  +  i/?(62+26-l).  In  the 
same  limit,  the  second  line  vanishes  (£=1+t— £,  so  that  £=1+t  when  C— 0)-  The  integral  occupying  the 
third  line  also  vanishes  then,  as  can  be  seen  by  dividing  the  interval  of  integration  (0,r)  into  (say)  the 
two  equal  parts  (0,|r)  and  (jT,t): 

|  [  T  dC  {[6  +  |ic(l  +  6)2]  Cos(kC^)  +  [1  +  |/c(l-«2)]  Sin(KC6)}[exp(-KC)]/^  I 
J  0  x 

<  [26+k(1+6)2]|  [exp(-K0]  dC/£  =  [26+«(l+6)2]|  |  [exp(-/<)]  -j-+  |  ^[exp(-/cC)]  -g- J 

<  [26+k(1+6)2]|  [exp(-KC)]d^  +  J  ^[exp(-KC)]  dC  j 

=  [26+/c(l+6)2]{exp(-Kr/2)  -exp(-fcr)  +  [l-exp(-Kx/2)]/(l+r/2)}//c  (A.8) 

The  last  line  of  Equation  (A.8)  vanishes  as  r- *oo;  so,  therefore,  must  the  integral  on  the  first  line 
(q.e.d.).  It  has  thus  been  shown  that  [Equation  (A.7)]: 

I(oo)  =  x[l  +  I/?(62+26-l)](a3/cp)P  =  AP  (A.9) 

Since  Poisson’s  ratio  v  fixes  the  constants  /?,  k  and  6,  Equation  (A.9)  dictates  that  if  two  systems  have 
values  of  v,  a  and  c^  that  make  A  the  same  for  both,  then  the  ratio  of  their  cavity- wall  pressures  Pj 

and  P2  will  equal  the  ratio  of  the  corresponding  total  impulses  I^oo)  and  I2(oo)  delivered  normal  to 

plane  £P  [I(oo)ng>  is  the  total-impulse  vector;  a op  =  unit  normal  to  <?].  That  is: 

I2(oo)/I1(oo)  =  P2/P1  (A.10) 

Systems  with  the  same  A- value  [so  that  Equation  (A.  10)  holds]  include  any  two  having  the  same  value 
of  both  v  and  a3/cp;  those  in  turn  include  any  two  with  the  same  values  of  v,  a  and  cp,  and  those 
include  any  two  with  the  same  values  of  G,  M,  a  and  material  density  p. 

On  that  same  space-time  domain  (i.e.,  the  domain  of  disturbance:  t  >  0,  0  <  r— a  <  cpt),  the 
radial  displacement  ur  is  related  to  radial  position  and  time  as  follows  (ibid.): 
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ur(£,r)  =  ur  =  { 1  -  [Cos(«C«)+|(l-20  Sin««)][  exp(-«C)]}  (A-n) 

Hence,  if  two  systems  have  the  same  values  of  G,  M  and  a,  then  the  ratio  of  their  displacements  at  any 
given  values  of  £  and  r  will  equal  the  ratio  of  their  cavity-wall  pressures;  indeed,  the  ratio  of  their 
(ur/a)-values  will  equal  the  pressure  ratio  if  G  and  M  alone  are  the  same  for  both.  Most  important,  if 
g'm,  a  and  p  are  the  same  in  the  two  systems,  then  [Equation  (A.10)]  the  total-impulse  ratio  is  also 

equal  to  the  ratio  of  pressures,  so  that 

K(£’r)]2/[ur(&r)]l  =  ^M/1^00)  (A'12^ 

Note  also  that,  with  cavity  radius  (a)  the  same  in  both  systems,  a  given  value  of  £  corresponds  to  the 
same  radius  r  for  both  (r=£a);  likewise,  with  M,  a  and  p  equal,  a  given  value  of  r  corresponds  to  the 

same  time  in  both  (t=ar  /cp;  Cp=M/ p). 
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APPENDIX  B 


CONTAINED  EXPLOSIONS  IN  SIMPLE  MATERIALS:  HOW  CAVITIES  FORM 

Myriad  conventional  charges,  and  most  nuclear  ones,  have  been  detonated  in  cavities  of  small  radius 
(0-.4  m/T1/3).  Displacements  of  near-cavity  material  are  then  appreciable  fractions  of  final  cavity 
radii.  The  theory  of  explosions’  effects  on  their  surroundings,  however,  is  better  developed  for  bursts  in 
so-called  “big  holes”  (>2  m/T1/3)  than  for  small  holes.  Least  well  understood  is  the  formation  of 
craters  by  bursts  shallow  enough  to  produce  them. 


In  brief,  following  explosive  energy  release  in  a  big  hole,  the  pressure  of  cavity  gases  can  still  be  large 
enough  to  bring  material  around  a  cavity  to  its  shear-strength  limits.  In  such  material,  further 
outward  displacement  of  the  cavity  wall  then  takes  place  without  further  increase  in  the  difference 
between  radial  stress  and  hoop  stress  (at  the  wall,  radial  stress  equals  cavity  pressure).  Since 
containers  (tires,  for  example)  hold  gas  under  pressure  only  by  dint  of  differences  between  radial  and 
hoop  stresses  in  their  walls,  part  of  the  pressure-load  that  would  be  opposed  by  hoop  stresses  in  elastic 
wall  material,  is  transferred  to  greater  distances,  where  outward  displacements  increase  as  a  result.  At 
some  medium-specific  distance  from  the  cavity,  and  beyond,  stresses  are  too  low  to  cause  strength 
limits  to  be  reached,  and  deformation  is  elastic.  As  a  result,  the  whole  system  comes  to  mechanical 
equilibrium  under  (a)  the  residual  pressure  of  cavity  gases  and  (b)  any  stresses  (e.g.,  overburden) 
present  throughout  undisturbed  material.  That  happens  sooner,  and  with  less  wall  displacement,  if  the 
hole  is  so  large  that  strength  limits  are  not  reached  at  all  in  adjacent  material. 

Patently  basic  to  big-hole  theory,  and  supported  by  its  experimental  validation,  is  the  following  idea. 
Geo-materials  deform  elastically  during  sufficiently  small  stress  excursions  about  overburden,  and  they 
flow  (i.e.,  move  with  shear  stresses  at  limiting  values)  during  larger  excursions.  For  bursts  in  small 
holes,  we  therefore  ask  what  results  follow  from  assuming  that  flow  at  high  stress,  and  elastic 
deformation  at  small  stress,  are  the  modes  of  behavior  that  dominate  cavity  growth.  The  relatively 
high  pressures  then  felt  by  surrounding  earth  cause  three  important  departures  from  the  motion  of 
material  around  a  big  hole:  (i)  An  appreciable  fraction  of  the  explosive  yield  is  converted  to  kinetic 
energy  of  the  earth,  (ii)  Displacement  of  the  cavity  wall  is  large  compared  to  its  initial  radius,  (iii) 
Flow  occurs  out  to  ranges  several  times  as  large  as  the  final  cavity  radius.  Quasistatic  and  small-strain 
approximations  suffice  to  describe  displacements  due  to  bursts  in  big  holes,  whereas,  in  view  of  (i)  and 
(ii),  the  full  equations  of  motion  for  arbitrary  strains  must  be  used  to  describe  cavity  growth  from 
bursts  in  small  holes.  Item  (iii)  holds  center  stage  here,  however:  Development  of  an  extended  region 
of  flow  profoundly  influences  not  only  the  size  of  the  final  cavity,  but  its  stability  and  the  state  of  the 
surrounding  medium,  as  will  now  be  described. 
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B.l  SIMPLE  MATERIALS. 

Even  when  performed  with  elaborate  material  models,  calculations  of  contained  bursts  in  small  holes 
almost  always  show  that  material  near  the  cavity  is  left  with  permanent  compressive  stresses  much 
larger  than  overburden.  These  stresses  have  their  origin  in  basic  properties  of  solids;  as  a  result,  their 
presence  can  be  understood  by  considering  simple  materials.  For  present  purposes,  we  therefore  label 
as  “simple”  any  isotropic  material  for  which,  in  a  spherical  field:  (A)  Volume  can  change  only 
elastically,  (B)  the  difference  between  radial  stress  oT  and  hoop  stress  ag  cannot  exceed  a  fixed  limiting 
value  Y  (its  “von  Mises  limit”),  and  (C)  deformation  is  elastic  as  long  as  the  resulting  values  of 
\<rr-<re\  remain  <  Y;  otherwise  the  condition  \<rT-<rg\='i  holds  during  deformation,  along  with 

condition  (A). 

As  background,  note  that  in  a  spherically  symmetric  field:  <rr  and  cg  are  principal  stresses;  shear  stress 
varies  with  the  orientation  of  the  surface  on  which  it  is  measured,  and  is  zero  on  coordinate  planes  of 
any  principal-axis  system;  the  shear  stress  at  a  point  is  proportional  to  aD  =  ar-ag,  and  its 
maximum  value  (||<rD|)  acts  on  any  plane  that  makes  a  45°  angle  with  the  radial  direction.  From 
these  facts,  it  follows  that  the  constraint  \<rD\  <Y  [condition  (C)]  does  no  more  (nor  less)  than  limit 
the  ability  of  simple  materials  to  support  shear  stress  -  in  accord  with  a  core  property  of  real  materials. 
Also,  if  the  shear  stress  attains  its  limiting  value  ^Y  in  a  simple  material,  it  does  so  on  planes  at  45  to 
the  radial.  Most  important:  If  the  shear  stress  on  such  a  plane  (its  maximum  at  any  given  point)  is 
maintained  at  its  limiting  value  in  a  simple  material,  then  any  slippage  that  occurs  along  the  plane 
generates  no  opposing  stresses.  By  contrast,  at  smaller  maximum  shear  stresses  than  that  limit,  even  a 
slight  additional  shearing  along  the  plane  will  cause  an  opposed  increase  in  shear  stress. 

B.2  RESIDUAL,  OR  “LOCKED-IN”  STRESS. 

Consider  now  the  behavior  of  simple  materials  in  uniaxial  strain  (which  spherical  deformation  becomes 
in  the  limit  of  large  radius).  In  particular,  with  <rr  and  ag  denoting  axial  and  lateral  stress, 
respectively,  let  an  element  of  simple  material  be  subjected  to  uniaxial  compression  (<rr>0)  and  then 
unloaded  back  to  <rr=0.  Loading  is  at  first  elastic  [condition  (C)];  <rg  will  increase  in  the  materials  of 
greatest  interest  (e.g.,  linear  elastic  materials  with  positive  Poisson’s  ratio,  v ),  but  not  as  much  as  <rT 
does  [except  in  fluids,  whose  shear  stress  is  always  zero  (at  least  quasistatically)].  Hence,  <rD 
(=<Tj.  _  <r g)  increases  with  <rr  until  its  limiting  value  Y  is  reached.  Further  growth  of  the  axial  load  <rT 
results  in  added  compression,  but  <rD  remains  equal  to  Y.  Thus,  once  <rD  reaches  the  value  Y, 
increases  in  compression  no  longer  give  rise  to  increases  in  shear  stress,  but  only  to  increases  in  mean 
stress  a  [  =  l(<rr+2ag)].  Yet,  further  axial  compression  is  automatically  attended  by  more  shear 
deformation.  In  contrast  with  earlier  loading  (when  <rD  <  Y),  however,  that  deformation  no  longer 
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generates  opposing  stresses.  Hence,  resistance  to  uniaxial  compression  drops  when  aD  attains  its 
limiting  value  Y;  a  compressive  increment  in  ar  causes  a  larger  decrease  in  volume  if  the  material  flows 
(|(T£j|=:Y)  than  if  it  deforms  elastically  (jcr^l  <  Y). 

Unloading,  which  is  initiated  when  ar  has  reached  a  maximum  value  £r,  relieves  shearing  stresses  (and 
others).  As  a  result  (since  shear  stress  is  proportional  to  cr^),  <rD  falls  from  its  limiting  value  Y  as 
unloading  begins,  and  unloading  proceeds  elastically.  As  already  noted,  however,  the  material  is  less 
compressible  elastically  than  when  it  flows,  Hence,  the  element  expands  less  when  <rr  is  reduced  from 
£r  to  zero,  than  it  shrank  on  increasing  crr  from  zero  to  £r.  The  net  effect  of  the  load-unload  cycle  is 
therefore  to  leave  the  element  with  residual  compression,  and  hence  [condition  (A)]  with  residual 
compressive  mean  stress.  Moreover,  since  <rr=0  after  unloading,  it  follows  that  the  final  value  of  <?$  is 
compressive.  Thus,  lateral  compressive  stress  remains  in  the  element  as  a  result  of  its  having  flowed 
during  compression.  Hereafter,  such  residual  stresses  are  termed  “locked-in  stresses. 

If  £  is  large  enough,  then  a  complication  to  the  load-unload  cycle  just  described  is  met.  \&j)\  W!U 
reach  its  limiting  value  on  unloading  Y).  In  particular,  just  as  in  elastic  loading,  <Tq  does  not 

change  as  rapidly  as  <rr  on  unloading,  so  that  \& p\  would  continually  increase  if  either  elastic  loading  or 
unloading  were  to  continue  indefinitely.  Once  the  limiting  shear  stress  -Y  is  reached  in  unloading, 
therefore,  <T£y=<Tr  —  aQ  must  remain  at  the  value  -Y.  Since  o-r=0  on  complete  unloading,  it  follows 
that  the  locked-in  stress  <Tq  cannot  exceed  Y.  For  materials  whose  elasticity  is  linear,  the  limiting 
locked-in  stress  Y  is  found  whenever  £r  >  <r *,  where  <=(cL/cs)2Y  and  the  material’s  longitudinal  and 
shear- wave  speeds,  respectively,  are  and  c^  [(cx/c5)2— M/G].  In  fact,  flow  then  starts  when  cr=^r^ 
and  the  locked-in  stress  increases  linearly  with  £r  from  zero  to  Y  for  values  of  £r  from  ^ <rr  to  <Jr.  For  a 
simple  material  whose  elasticity  is  linear,  Figure  B-l  illustrates  the  stages  of  a  load-unload  cycle  of 
uniaxial  strain. 


As  for  spherical  fields,  hoop  strains  play  a  major  role  in  their  development  (B.3  below).  As  a  result, 
stress  paths  depart  somewhat  from  those  of  uniaxial  strain  during  compressive  loading,  and  on 
subsequent  unloading  the  differences  are  large  (B.4).  Still,  when  a  simple  material  is  driven  by  a 
shock  or  a  steeply  rising  compression  pulse  to  states  in  which  a  >>Y,  a  locked-in  hoop  stress  about 
equal  to  Y  can  be  expected  after  unloading.  Accordingly,  to  the  extent  that  elasticity  and  flow  exhaust 
the  modes  of  deformation  of  material  near  an  actual  explosively  formed  cavity,  that  material  should 
afterward  contain  locked-in  stresses  of  a  size  comparable  to  its  strength.  For  such  events, 
computational  models  of  geo-materials  are  not  “simple,”  but  elasticity  and  flow  do  dominate  their 
calculated  behavior  near  the  cavity  -  and  the  models  do  yield  radial  distributions  of  locked-in  stress  of 
the  expected  size.  We  now  recall  how  those  distributions  evolve. 
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B.3  THE  “ATLAS  ZONE”  AND  ITS  DEVELOPMENT. 

As  calculated  using  diverse  material  models  (many  quite  complex)  and  discretized  equations  of  motion 
(but  otherwise  exact),  cavity  growth  ends  for  contained  bursts  at  about  2-  (but  less  than  4)  ms/T  ^  . 
The  final  cavity  radius  greatly  exceeds  the  value  it  would  attain  by  quasistatic  expansion  (B.5,  below) 
and  the  final  cavity  pressure  is  often  smaller  than  overburden.  Hence,  the  cavity  is  not  held  open  by 
the  pressure  of  gases  within.  Instead,  the  forces  that  hold  it  open  in  stable  equilibrium  with  the 
surrounding  earth  are  supplied  by  locked-in  stresses.  Those  stresses  develop  during  cavity  growth  over 
the  region  of  inelastic  deformation,  which  covers  a  spherical  shell  of  finite  thickness  bounded  by  the 
cavity  wall.  Since  its  stresses  support  the  load  of  an  adjoining  world,  that  shell  (more  precisely,  an 
inner  core  thereof;  B.5)  is  referred  to  here  as  the  “Atlas  Zone”  To  see  how  an  Atlas  Zone  forms, 
consider  a  tamped,  spherically  symmetric  explosion  in  homogeneous  simple  material,  under  a  uniform 
overburden  great  enough  to  preclude  cracking.  Development  of  an  Atlas  Zone  then  takes  place  as 
follows. 

Pressures  of  10-to-20  GPa  from  conventional  explosives  (a  hundred  times  that  in  nuclear  bursts)  drive 
nearby  material  away  from  the  burst  point  at  high  speed,  as  cavity  growth  begins.  As  a  result,  when 
the  pressure  in  the  cavity  falls  to  overburden  levels,  material  near  the  wall  is  still  moving  outward. 
Further  deceleration  of  wall  material  then  occurs,  as  before,  because  cavity  pressure  continues  to  drop, 
and  because  radial  stresses  are  more  compressive  than  hoop  stresses  [which  makes  for  negative 
acceleration;  Equation  (2.1)]. 

The  stress  difference  <rD  (  =  crr  -  ae)  remains  at  its  von  Mises  limit  Y  through  almost  the  entire  time 
of  cavity  growth,  despite  the  drop  in  cavity  pressure,  owing  to  continually  increasing  hoop  strains. 
That  is,  expansion  of  the  cavity  causes  rings  of  material  around  it  to  stretch,  constantly  driving  hoop 
stresses  toward  tensile  values,  whereas  radial  stresses  at  the  cavity  wall  must  equal  the  (compressive) 
cavity  pressure.  The  difference  & jj  between  radial  and  hoop  stresses  is  then  constant  at  Y>0,  and 
hence  acts  [Equation  (2.1)]  as  a  source  of  inward  acceleration—  a  constant  deceleration  as  material 
slows  (ignoring  small  fractional  changes  in  its  radius).  Once  velocity  becomes  inward-directed, 
however,  hoop  strains  and  stresses  grow  more  compressive  and  <t d  decreases  from  Y:  As  inward 
motion  proceeds,  material  deforms  elastically  near  the  cavity  wall. 

Now,  in  states  of  shear  failure,  unlimited  deformation  can  occur  with  no  change  in  stress  (B.l),  but 
elastic  deformation  is  always  attended  by  opposing  stress-changes.  Inward  motion  therefore  gives  rise 
to  ever-growing  stresses  to  oppose  that  motion  (in  salt,  for  example,  an  elastic  decrease  in  cavity  radius 
by  1%  would  produce  a  compressive  hoop-stress  increase  of  a  few  tenths  of  a  GPa).  As  a  result,  when 
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the  straining  that  causes  inward  acceleration  of  wall  material  finally  arrests  its  outward  motion,  no 
more  than  a  slight  decrease  in  cavity  radius  can  occur  without  generating  compressive  hoop  stresses 
that  would  impart  large  outward  acceleration  to  the  wall.  In  the  event,  the  wall’s  velocity  falls 
continuously  to  zero  and  then  turns  negative.  The  resulting  compressive  hoop-stress  increments 
-  locked-in  stresses  of  the  kind  discussed  above  (B.2)  -  therefore  simply  arrest  its  inward  motion  before 
any  large  inward  velocity  develops:  After  reaching  its  maximum  radius,  the  wall  oscillates  with 
negligible  (and  decreasing)  amplitude  about  its  final  equilibrium  position. 

B.4  LOCKED-IN  STRESS:  THE  SALMON  EVENT. 

As  already  noted  (B.2,  above),  the  stress  paths  traced  in  spherical  bursts  differ  substantially  from  the 
paths  (Figure  B-l)  followed  in  uniaxial  strain.  In  particular,  for  reasons  stated  in  Section  B.3,  the  wall 
and  nearby  earth  continue  to  flow  as  they  slow  down.  Hence,  as  pressure  falls  in  the  expanding  cavity, 
unloading  occurs  largely  along  the  limiting  stress  path  aD  =  cr  -  <re=Y:  Loading  stress-paths  (which 
do  resemble  that  of  Figure  B-l  for  uniaxial  strain)  are  in  large  part  retraced  on  unloading  (strain  paths 
are  not).  With  the  onset  of  inward  wall  motion  (negative  radial  velocity),  further  unloading  takes  place 
elastically  and  with  almost  constant  pressure  at  the  cavity  wall  (B.3).  Stress  paths  for  elastic 
unloading  at  constant  crr>0  are  therefore  parallel  to  the  line  <rr=0  in  Figure  B-l  —  and  final  states  are 
reached  by  material  near  the  wall  along  paths  parallel  to  that  one. 

Since  the  wall  reverses  direction  sooner  than  material  further  out,  an  elastic  unloading  wave  spreads 
outward  through  the  region  of  flowing  material.  The  development  of  large  compressive  hoop  stresses 
requires  little  inward  motion,  however  (B.3),  and  stops  such  motion  after  only  slight  inward 
displacement.  Radial  stresses  are  larger  in  the  interior  of  the  flow  region  than  at  the  wall  (B.5;  also, 
falling  gas  pressure  is  felt  first  at  the  wall),  and  hence  are  not  fully  relieved  by  that  displacement. 
Thus,  as  motion  winds  down  near  the  cavity,  a  characteristic  stress-pattern  emerges,  of  which  an 
example  is  discussed  in  the  rest  of  this  section:  Hoop  stress  has  its  maximum  compressive  value  at  or 
very  near  the  cavity  wall,  and  then  decays  smoothly  with  radius  through  the  flow  region.  Radial 
stress,  which  is  equal  at  the  wall  to  the  final  gas  pressure,  rises  to  a  maximum  at  a  somewhat  larger 
radius;  spherical  divergence  (presumably)  causes  stress-decay  at  still  larger  radii  (geometric  attenuation 
of  displacement,  stress,  etc.,  with  radius  is  especially  rapid  on  the  region  of  flow,  which  abuts  the 
cavity  and  contains  all  locked-in  stress). 

The  processes  just  described,  and  also  those  of  Section  B.3,  are  evident  in  stress  and  radial-velocity 
wavetrains  (Figures  B-2  and  B-3,  respectively)  from  an  ATI  calculation  of  the  Salmon  event.  As 
modeled  therein,  the  strength  of  dome-salt  increases  smoothly  with  mean  stress  W  (from  4  MPa  at  a  =0 
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Figure  B-2.  Stress  vs.  slant  range,  calculated  for  the  Salmon  Event  37.0,  45.4  and  70.0  milliseconds 
(ms)  after  the  burst.  Radial  and  hoop  stresses  are  shown,  respectively,  by  curves  marked 
“x”  and  “0”;  unmarked  curves  refer  to  mean  stress. 
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to  41  MPa  at  W  =50  MPa,  to  87  MPa  for  a  >  220  MPa);  the  model  also  admits  inelastic  changes  in 
volume.  In  the  calculation,  the  cavity  wall  moves  outward  for  ~  42  ms  after  its  motion  starts.  At  37 
ms,  the  medium  is  flowing  at  all  points  between  the  wall  (Rc  «  16  m  and  moving  outward  at  ~  25 
m/s;  Figure  B-3)  and  a  local  peak  almost  nine  cavity  radii  further  out  (near  the  wavefront;  Figure  B- 
2).  By  45.4  ms,  the  wall  is  moving  inward  (at  ~  7  m/s),  and  material  is  unloading  elastically  from  Rc 
(the  wall-radius)  out  to  ~  2RC.  At  70  ms,  near-static  equilibrium  has  been  reached  from  Rc  (still  «  16 
m)  to  6|RC.  Evidence  of  that  lies  in  small-to-negligible  radial  velocities  over  the  region  noted.  More 
compelling:  At  static  equilibrium,  <rr  must  equal  <re  wherever  a<rr/ar=0  [Equation  (2.1)],  and,  in 
Figure  B-2  at  70  ms,  hoop  stress  does  equal  radial  stress  at  a  point  very  close  to  the  latter’s  maximum 
(defined  here  as  the  outer  radius  of  the  Atlas  zone;  Section  B.3).  Thus,  the  figure’s  peak  hoop  stress 
(  ~  .100  GPa  above  an  overburden  of  18  MPa)  is  locked  into  the  medium. 

For  particles  of  salt  that  end  up  at  §RC,  4RC  and  8RC,  Figure  B-4  presents  the  stress  paths  traced  in 
reaching  final  states  of  stress.  Evidently,  after  a  brief  elastic  stage,  loading  takes  place  at  the  limit  of 
shear  stress  (aD  =aT-ae= Y),  which,  on  this  model,  is  reduced  by  shock  heating  (such  reduction,  a 
slight  one,  can  be  seen  on  path  §RC,  but  not  for  more  distant  material).  Further,  until  the  final  stress- 
state  is  approached,  most  unloading  also  takes  place  mainly  at  that  limit;  path  4RC  departs  from  it 
somewhat  (as  do  other  paths)  when,  for  a  short  time  after  the  initial  rise  to  peak  compression 
(loading),  cT  decreases  faster  than  (r0.  Soon  after  material  starts  moving  inward,  however,  its 
unloading  proceeds  elastically.  Near  the  cavity  wall,  where  <rr  is  about  equal  to  cavity  pressure  P,  that 
happens  at  once:  With  P  (hence  <rT)  much  less  sensitive  than  hoop  stress  <r0  to  elastic  decreases  in 
wall-radius,  er D  falls  from  its  shear-failure  value  Y  if  elasticity  is  assumed  [even  though  inward  motion 
of  the  wall  makes  both  cr6  and  P  (hence  <rr)  less  tensile]  -  whence  changes  in  strain  must  be  elastic. 
Elastic  behavior  ensures  a  large  drop  in  <r jj  without  much  inward  displacement  —  hence  little  change  in 
P  -  and,  perforce,  in  <rr.  Elastic  unloading  (its  final  stage)  therefore  unfolds  near  the  cavity-wall  in  the 
general  direction  of  lines  of  constant  <rr  (Figure  B-l),  as  Figure  B-3  affirms. 

B.5  CAVITY  SIZE  AND  ELASTIC  TRANSPORT:  DYNAMIC  AND  GEOMETRIC  EFFECTS. 

The  fact  that  large  strains  develop  near  contained  spherical  bursts  in  small  holes  (0-.5  m/T1/3  radius), 
follows  from  the  cavities’  large  ratios  of  final  radius  to  initial  radius.  By  the  same  token,  strains  are 
large  even  around  quasistatically  pressurized  cavities  (to  burst-like  pressures).  Hence,  cavity  growth  is 
inherently  nonlinear,  as  shown  by  the  quasistatic  elasto-plastic  result  that  final  cavity  radius  grows 
exponentially  with  cavity  pressure  (above  the  least  pressure  at  which  the  wall  fails  in  shear;  Hill,  loc. 
cit.).  Quasistatic  release  of  explosive  energy  in  a  small  cavity  results  in  a  final  radius  of  ~  1.19  Y~ 
m/T1/3,  where  Y,  in  MPa,  is  the  maximum  of  Y(o )  [e.g.,  87  MPa  for  the  salt  model  of  Section  B.4; 
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Figure  B-4.  Calculation  of  the  Salmon  Event:  Stress-paths  followed  by  sample  particles  at  final  radii 
equal  to  1.5,  4  and  8  times  the  cavity  radius,  Rc.  In  all  cases,  compressive  loading  starts 
from  a  hydrostatic  state  at  .181  kb  of  mean-stress  (the  overburden),  and  quickly  takes  the 
shear-stress  (and  hence  the  stress  difference  <rD)  to  limiting  values,  where  they  remain  until 
loading  ends.  At  1.5  Rc  and  4  Rc,  peak  values  of  mean  stress  are  to  large  to  be  shown  on 
the  plot.  Most  subsequent  unloading  takes  place  with  u jj  retracing  the  same  limiting 
values  traversed  on  loading.  In  its  final  stage,  unloading  is  elastic,  and  radial  stresses 
increase  at  1.5  Rc  and  8  Rc.  In  the  computational  salt  medium,  shear  strength  (and 
hence  limiting  values  of  <td)  vary  with  internal  energy  as  well  as  mean  stress;  thermal 
effects  on  strength  are  insignificant  at  4  Rc  and  8  Rc,  but  noticeable  at  1.5  Rc.  Also, 
inelastic  volume  increases  (“bulking”),  attend  shear  failure,  while  inelastic  decreases  in 
volume  (“compaction”)  take  place  in  hydrostatic  compression. 
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writing  Y(F)  distinguishes  strength-as-a-function-of-mean-stress-er  from  its  maximum  value  Y]. 
Dynamic  effects  (due  to  explosive  release  of  energy)  increase  the  final  radius  to  ~  2.86  Y"2/7  m/T1/3. 

The  first  of  these  two  expressions  fits  closely  an  accurate  pencil-and-paper  solution  to  the  quasistatic 
problem.  The  second  is  a  close  fit  to  results  obtained  by  massive  computing  (using  discretized 
equations  of  motion)  for  spherical  bursts  in  simple  materials.  So  far,  it  has  nonetheless  proven  accurate 
to  within  ±  25%  for  media  with  porosities  <  5%  and  Y  <  170  MPa,  with  above-average  radii  obtained 
(usually)  in  the  fields  computed  for  non-simple  materials.  More  important,  final  cavity  radii  measured 
for  contained  tamped  explosions  are  also  larger -about  1.0  m/T1/3  for  ordinary  and  nuclear  explosives 
(the  formula  gives  .80  m/T1/3  for  dome-salt).  At  that  radius,  accepted  equations  of  state  for  cavity 
gases  (vaporized  earth,  in  essence,  for  nuclear  bursts)  yield  pressures  of  ~  .28  GPa  if  all  energy  stays  in 
the  cavity;  massive  computing  leaves  |  to  |  of  the  energy  in  the  cavity  when  it  stabilizes,  which 
implies  that  final  cavity  pressures  are  generally  <  100  MPa  (chemical  explosive)  and  50  MPa  (nuclear). 

Seconds,  it  should  be  noted,  are  long  compared  to  the  times  called  “final”  here  (Section  B.4);  over 
much  longer  times  still,  cavity  gases  cool  (or  may  diffuse  into  nearby  soil/rock)  and  their  pressures 
drop  by  large  factors.  Salmon  presents  a  dramatic  case  in  point:  Re-entry  drilling  disclosed  a  cavity 
whose  pressure,  six  months  after  the  event,  had  fallen  to  .03  MPa  —  a  partial  vacuum  inside  a  cavity  of 
16.7-m  radius  (.96  m/T1/3;  Rawson,  et  al.,  1966).  Moreover,  at  times  of  a  second  or  less,  an 
appreciably  larger  radius  (by  a  factor  of  almost  1.5)  has  been  postulated  for  the  Salmon  cavity 
(reduced  by  creep  to  16.7  m  over  six  months’  time;  Rogers,  1966).  As  estimates  of  actual  “final” 
pressures,  values  up  to  100  MPa  (chemical)  and  50  MPa  (nuclear)  could  therefore  be  excessive.  On  the 
other  hand,  <rr  attains  values  much  greater  than  <r*  near  the  wall  (Section  B.3),  because  wall-pressure 
does -and,  in  media  well  described  as  simple,  locked-in  stresses  therefore  reach  the  maximum  value 
possible  in  the  initial  shock  (i.e.,  in  uniaxial  strain).  That  value,  Y,  most  likely  exceeds  the  “final 
cavity  pressure.  Whether  it  does  or  not,  however,  geometry  magnifies  the  part  played  by  locked-in 
stress  (relative  to  gas  pressure)  in  holding  open  a  final  spherical  cavity  -  as  we  now  explain. 

The  instructive  geometric  object  (Figure  B-5)  is  part  of  a  right  circular  cone  of  small  angle  6  with  its 
vertex  at  the  burst  point  -  a  frustum,  in  the  limit  6-*0.  For  finite  9,  the  cone-segment  'S  in  question 
runs  from  the  (spherical)  cavity-wall  at  its  “final”  radius  Rc  to  the  spherical  sheet  of  material  whose 
radius  Rp  at  “final”  equilibrium  bounds  the  region  of  flow.  That  is,  end-caps  —  sections  of  spherical 
surfaces  centered  at  the  cone’s  vertex  -  truncate  the  cone  to  form  IP;  further,  material  at  Ts  outer  end- 
cap  has  undergone  shear  failure  and  ?  holds  all  cone  material  that  ever  has  so  failed.  Symmetry  makes 
the  normal  to  Ts  surface  at  any  point  a  principal-stress  (and  -strain)  axis,  so  that  traction  thereon  is 
normally  directed.  In  the  direction  of  the  cone’s  axis,  which  is  radial,  the  force  on  an  areal  element  of 
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Ts  surface  is  thus  equal  to  the  element’s  area,  times  the  normal  traction  there  (hoop  stress  over  the 
conical  surface,  radial  stress  on  the  end-caps),  times  the  axial  (i.e.,  radial)  component  of  the  inner 
normal  to  the  element.  In  turn,  the  product  of  the  normal’s  axial  component  with  the  element’s  area  is 
just  the  area  of  the  element’s  projection  on  any  plane  normal  to  the  cone’s  axis. 

The  projection  of  the  inner  end-cap  on  such  a  plane  is  clearly  a  circle  of  radius  RcSin|  [and  area 
TrR^(Sinf)2];  since,  over  the  whole  cap,  normal  traction  equals  final  wall  pressure  P ,  the  radial  force 
Fc  on  that  cap  is  P'7rR^(Sin|)2.  In  the  same  plane,  the  projection  of  if’s  conical  surface  is  evidently  an 
annulus  of  inner  radius  RcSinf  and  outer  radius  RFSinf,  but  hoop  stress  (normal  traction)  thereon 
varies  with  radius;  accordingly,  the  radial  force  FL  delivered  to  the  conical  surface  is  the  product  of  its 
projected  area  jr(R^-R^)(Sinf)2  with  an  average  (?',  say)  of  the  values  of  final  hoop  stress  over 
said  surface.  Thus,  as  the  ratio  (<r^/  P  ) ( RF - K2. )  /  R  of  the  two  radial  forces  shows,  the  geometric 
factor  (Rf-R£.)  /R^  magnifies  the  contribution  of  locked-in  hoop  stress  to  the  outward  force  on  the 
cone,  relative  to  that  of  cavity  pressure;  the  factor  is  3  if  RF  is  just  twice  Rf.  but  it  can  be  much 
greater  [for  Salmon  (Figure  B-3,  bottom)  actually  exceeds  P'  out  to  ~  3RC,  and  RF  «  10  Rc]  . 

The  force  Ff  (radially  inward)  exerted  on  the  outer  end-cap  is  equal  to  -<rFirRF(Sin|)2,  where  <xE 
denotes  the  final  radial  stress  at  RF  (the  elastic-plastic  radius);  FL,  expressed  exactly  in  terms  of  <r~, 
equals  d(rSinf)2=7r(Sinf)2/^dr2,  over  Rc  <  r  <  RF.  Since  Fc+Fl+Ff  must  equal  zero  at 

equilibrium,  the  force  -Ff  [=<r^xR|.(Sinf)2]  delivered  to  unfailed  cone-material  across  the  elastic- 
plastic  boundary  (!f’s  outer  end-cap,  at  radius  RF)  must  equal  F^.+Ff.  Given  purely  elastic 
deformation  beyond  RF,  final  radial  displacement  there  is  determined  by  o'E  (plus  the  overburden  at 
huge  radii  where,  physically,  the  burst  cannot  effect  finite  changes)  —  and,  <rF=  (Fc+Fl)/R^.  Indeed, 
the  displacement  at  any  r  >  lip  is  proportional  to  <x ^  if  the  state  of  zero  stress  is  defined  as  that  of  the 
pre-shot  medium  (whereupon  overburden  is  zero  by  definition);  cr~  (for  example)  is  then  automatically 
the  permanent  change  in  radial  stress,  at  a  given  field-point,  wrought  by  the  burst  [initial  stresses  are 
otherwise  subtracted  from  cr etc.,  in  the  formulas  above  -  which  still  remain  valid  because,  over  a 
closed  surface  in  any  equilibrium  field  (including  the  initial  field),  the  net  force  is  zero]. 

The  spherical  sheet  of  material  that  ends  up  at  Rp  could,  mathematically  at  least,  have  started  out  at 
a  radius  >  Rp.  Given  purely  elastic  deformation  outside  that  sheet,  however,  the  change  in  radial 
stress  at  Rp  would  then  be  zero  or  tensile,  even  though  the  volume  of  solid  between  the  cavity  and  Rp 
has  decreased  and  all  volume  changes  are  elastic.  On  its  face,  such  an  outcome  flouts  physical  sense, 
which  lies  instead  (for  example)  with  an  exact  quasistatic  result  for  simple  materials:  At  the  elastic- 
plastic  boundary,  the  change  in  radial  stress  (|Y;  Hill,  loc.  cit.)  is  compressive.  As  for  Salmon,  near- 
equilibrium  obtains  at  the  last  time  shown  (Figure  B-3)  only  out  to  r«  6R^,  whereas  Rp  «  10  R^. 
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Still,  to  61^,  we  find  <r’> 0  (i.e.,  hoop-stress  changes  are  compressive),  and  the  same  holds  for  radial 
stress  (Figure  B-2);  results  of  the  calculation  at  later  times  (not  shown)  have  <t*>0  out  to  RF.  Thus, 
not  only  is  FL>0  (indeed,  jV' dr2  increases  steadily  with  r),  but  geometry  renders  Fc  (also  >0)  «  FL. 
The  change  <r£  in  radial  stress  at  RF  must  again  be  compressive  (since  Fc,+Fl>0),  so  that  material 
comes  to  rest  there  after  a  net  outward  displacement;  the  computed  field  directly  confirms  these  things. 

The  great  resistance  of  simple  media  to  elastic  changes  in  hoop  strain  provides  general  reason  to  expect 
such  results.  True,  inward  motion  at  and  near  the  cavity  wall  is  opposed  by  increasing  pressure 
thereon;  the  modulus  37P  (3P  to  4P)  relates  increments  in  wall-pressure  P  (hence  radial  stress)  and 
hoop-strain  [Cooper,  1966,  Equation  6].  By  contrast,  with  <r=  P,  the  modulus  relating  increments  in 
hoop-stress  and  -strain  is  37PA/M-KM,  where  £  is  >  1  for  0  <  v  <  i  [e  =  (l-2i/)(l+i')/(l-*')2)*  With 
P<.1  GPa,  M  far  exceeds  37P  [Salmon:  M/(37P)>100].  Hence,  inward  motion  (elastic  unloading)  near 
the  wall  is  resisted  mainly  by  hoop-stress  growth;  locked-in  hoop  stress,  not  gas  pressure,  keeps  the 
cavity  open.  Further,  as  inward  motion  is  arrested  and  equilibrium  impends,  that  sharp  growth 
(relative  to  radial  stress)  takes  hoop  stress  to  compressive  levels  above  those  of  radial  stress  (up  to  the 
limiting  stress-difference  Y).  At  the  same  time,  spherical  convergence  during  inward  motion  drives  up 
mean-stress  values,  which  tends  to  increase  radial  stress  near  the  cavity  -  except  at  the  wall,  where  <rT 
must  equal  P.  The  outlines  of  an  equilibrium  stress-distribution  on  the  flow  region  thus  emerge:  Hoop 
stress  rises  to  a  peak,  and  exceeds  radial  stress,  very  near  the  wall;  radial  stress  increases  more 
gradually  with  radius,  but  must  eventually  exhibit  the  decay  that  spherical  divergence  demands -and, 
at  its  peak,  ar  must  equal  <Je  (Section  B.4). 

At  the  opposite  extreme  (big  holes,  full  decoupling),  no  flow  takes  place  at  all  (RF=RC).  Particle 
displacements,  including  that  of  the  cavity  wall,  then  result  solely  from  a  permanent  change  in  cavity 
pressure.  In  all  cases,  though,  Fc  (gas  pressure)  and  (locked-in  stress)  determine  radial  stress,  and 
hence  displacement,  at  the  elastic-plastic  boundary  (r=RF).  In  turn,  that  displacement  determines 
how  much  of  the  cavity’s  volume  is  transported  elastically  to  distant  parts -just  as  mean  stress  locked 
into  the  flow-region  fixes  the  volume  added  to  the  cavity  by  elastic  compression  of  solid  on  that  region. 
The  ratio  Ff/(F(7+Ff)  measures  the  importance  of  locked-in  stress  to  the  final  field  of  displacement. 
Its  value  is  w  1  for  tamped  bursts  [>1  if  P’<0  (final  gas  pressure  <  overburden)  -  no  miracle],  and 
decreases  to  zero  as  hole-size  increases  to  the  minimum  level  for  full  decoupling. 

In  sum:  Residual  gas  pressure  dictates  the  displacement  field  left  by  an  explosion  in  a  large  hole.  For 
bursts  in  small  holes,  residual  (locked-in)  stress  in  the  host  medium  determines  permanent 
displacement.  Calculations  set  the  radii  of  “small”  holes  at  0-to-|  m/T1/3;  holes  with  radii  >2  m/T  ^ 
are  “large.”  These  results  hold  if  inelastic  volume  decreases  on  the  region  of  flow  aggregate  to  a  small 


B-14 


fraction  of  final  cavity  volume,  and  if  deformation  is  elastic  outside  the  region  of  flow.  That  is,  they 
hold  for  the  geologic  media  envisioned  by  big- hole  theory. 
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